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ABSTRACT 


The  infrared  radiance  arising  from  the  earth's  limb  for  altitudes  from  sea 
level  to  500  kilometers  and  within  the  spectral  limits  of  5  to  25  microns  was 
determined  from  theoretical  calculations.  The  atmosphere  was  divided  into 
two  altitude  regimes:  the  lower  atmosphere  (below  70  kilometers )  and  the 
upper  atmosphere  (above  70  kilometers).  The  lower  limb  models  were 
developed  from  the  University  of  Michigan  computer  program  with  modifica¬ 
tions,  including  the  addition  of  a  model  for  nitric  acid.  Spectral  radiances 
were  determined  as  a  function  of  geometrical,  geographical,  and  temporal 
parameters.  Emphasis  was  placed  on  developing  limb  radiance  models  for 
the  upper  atmosphere,  including  identification  of  active  molecular  species  in 
the  infrared,  the  specification  of  their  abundances  and  their  wavelengths  of 
emission.  In  addition,  particulate  material  sources,  including  infrared  radia¬ 
tion  arising  from  dust  particles  and  noctilucent  clouds,  and  uncertainties  in 
atmospheric  parameters  were  considered. 


in 
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FOREWORD 


This  special  scientific  repoi't  presents  the  results  of  a  theoretical  stud}  of 
the  emission  properties  of  the  earth's  atmosphere.  The  study  was  performed 
jointly  under  contract  F19628-69-C-0258  by  Honeywell  Inc.  as  prime  contractor 
and  GCA  Corporation  as  subcontractor. 

The  study  investigated  the  infrared  radiance  arising  from  the  earth's  limb 
from  sea  level  to  an  altitude  of  500  kilometers  and  within  the  electromagnetic 
spectral  limits  of  5  to  2  5  microns.  The  primary  purpose  of  the  study  was  to 
obtain  first-order  numerical  estimates  of  horizon  radiance  for  the  design  of 
a  horizon  experiment  program.  In  many  cases,  the  results  were  obtained 
without  experimental  verification;  users  are  therefore  cautioned  on  the 
limitations  of  the  resulting  data.  For  study  purposes,  the  atmosphere  was 
li\  idl'd  into  two  altitude  regimes:  the  lower  atmosphere  (below  70  kilometers) 
and  the  upper  atmosphere  (above  70  kilometers). 

Th  >  lower  limb  models  were  developed  from  the  University  of  Michigan 
computer  program  with  modifications,  including  the  addition  of  a  model  for 
nitric  acid.  Spectral  radiances  were  determined  as  a  function  of  geometrical, 
geographical,  and  temporal  parameters,  taking  into  account  the  range  of 
carbon  dioxide,  water  vapor  and  ozone  concentration,  and  temperature  profile 

variations . 

Emphasis  was  placed  on  developing  limb  radiance  models  for  the  upper  atmos¬ 
phere.  This  portion  of  the  study  included  the  identification  of  active  molecular 
species  in  the  infrared,  the  specification  of  their  abundances,  their  wavelengths 
of  emission  and  the  development  of  a  realistic  physical  model  to  perfor  m 
theoretical  radiance  calculations.  In  addition,  particulate  material  sources, 
including  infrared  radiation  arising  from  dust  particles  and  noctilucent 
clouds,  and  the  variability  of  limb  radiance  due  to  the  means,  extremes,  and 
uncertainties  in  atmospheric  parameters  were  considered. 
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SECTION  I 
INTRODUCTION 


This  study  is  concerned  with  the  problem  of  investigating  the  natural  infrared 
background  of  the  earth  in  the  5  to  25  micrometer  region  of  the  electromag¬ 
netic  spectrum.  This  general  problem  area  is  of  interest  from  the  standpoints 
of  systems  design,  military  surveillance  and  the  advancement  of  scientific 
knowledge.  The  estimates  of  the  infrared  radiant  background  from  physical 
models  developec  in  the  study  are  useful  in  these  areas.  Knowledge  of  the 
magnitude,  extreme  range  and  variability  of  the  background  radiances  is 
important  in  defining  the  specifications  of  radiometer  systems  and  the  logic 
of  developing  optimum  background  measurement  programs.  In  a  similar 
manner,  the  background  estimates  are  of  interest  in  the  areas  of  military 
surveillance  and  homing  interception  defense  measures.  F  inally,  the  problem 
area  is  extremely  interesting  in  terms  of  our  scientific  knowledge  of  the  physi¬ 
cal  processes  taking  place  in  the  atmosphere,  especially  in  the  altitude  region 
above  about  70  km. 

The  scope  of  the  current  effort  is  the  development  of  physical  models  to  esti¬ 
mate  infrared  radiation  arising  from  active  molecular  gases  and  particulate 
material  in  the  atmosphere  from  the  surface  to  about  500  km.  The  effort  is 
confined  to  viewing  the  atmosphere,  geometrically  speaking,  against  the  rela¬ 
tively  "cold"  background  of  space,  and  viewing  geometries  include  exoatmos- 
pheric  viewing  where  the  reference  frame  is  outside  the  earth's  atmosphere 
(e.g.,  viewing  the  atmospheric  limb  from  a  satellite  platform)  and  endoatmos- 
pheric  viewing  where  the  reference  frame  is  imbedded  within  the  atmosphere 
(e.g.,  viewing  from  the  zenith  to  the  horizon  from  a  rocket  platform).  This 
study  does  not  consider  the  infrared  background  of  the  earth-atmosphere 
system,  where  the  atmosphere  is  viewed  against  the  "hot"  infrared  source  of 
the  earth's  surface,  since  the  problem  area  has  been  investigated  extensively 
by  other  researchers. 
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The  basic  physics  of  developing  a  realistic  theoretical  model  to  estimate 
radiance  in  the  infrared  spectrum  are  reasonably  well  known.  Any  theoretical 
model  relies  on  a  knowledge  of  atmospheric  properties.  Our  knowledge  of  the 
physical  properties  of  the  atmosphere  such  as  temperature  and  composition 
are  quite  certain  from  the  surface  to  upper  balloon  altitudes  of  about  30  km, 
are  less  certain  in  the  low  altitude  rocket  range  from  30  to  70  km,  and  become 
uncertain  and  sometimes  speculative  at  altitudes  about  70  km. 

In  much  the  same  manner,  the  physical  methods  of  treating  the  problem  of 
computing  background  radiances  may  be  divided  into  two  altitude  regions,  a 
lower  atmosphere  from  the  surface  to  70  km,  and  an  upper  atmosphere 
above  70  km.  An  extensive  amount  of  theoretical  modeling  effort  has  been 
performed  in  the  lower  atmospheric  region.  Much  less  effort  has  been  per¬ 
formed  on  the  upper  atmospheric  region.  The  primary  effort  of  the  current 
study  is  directed  toward  the  problem  of  developing  a  realistic  physical 
model  in  the  upper  atmosphere. 

For  convenience  in  subsequent  discussions,  the  atmosphere  is  divided  into 
two  regions,  the  lower  atmosphere  and  the  upper  atmosphere.  Some  of  the 
physical  properties  distinguishing  these  two  regions  are  as  follows: 

Lower  Atmosphere  (surface  to  70  km) 

•  Physical  properties  such  as  chemical  composition,  temperature 
and  pressure  are  reasonably  well  known. 

•  Collision  between  molecules  is  rapid  enough  to  maintain  a  thermal 
equilibrium  population  of  vibrational  states. 

•  Radiative  transfer  theory  is  well  known  and  has  been  verified 
against  observations. 
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•  A  mixed  Lorentz-Doppler  line  shape  is  required  in  the  treatment 
of  atmospheric  transmission  properties. 

Upper  Atmosphere  (above  70  km) 


•  Atmospheric  properties  such  as  chemical  abundances  and 
temperature  are  uncertain  as  to  their  magnitude,  variability 
and  even  as  to  their  identity  in  the  case  of  minor  atmospheric 
constituents. 

•  Collisions  between  molecules  are  less  effective  in  maintaining 
excited  vibrational  states  as  compared  to  the  lower  atmosphere. 

•  Radiation  excitation  mechanisms  become  more  important  with 
increasing  altitude. 

•  Radiative  transfer  theory  is  reasonably  well  known  but  has  not 
been  verified  extensively  against  observations. 

•  A  Doppler  line  shape  is  a  good  approximation  in  treating  atmospheric 
transmission  properties. 

•  Chemiluminescent  and  fluorescent  processes  become  important. 

In  the  problem  of  investigating  the  natural  infrared  background  of  the  earth, 
the  current  study  relies  on  the  physical  model  developed  by  Anding  (1968)  in 
order  to  perform  radiance  calculations  for  the  lower  atmosphere.  This  model 
has  been  incorporated  in  a  working  computer  program  whose  calculations  have 
been  verified  against  observations.  The  computational  model  for  the  upper 
atmosphere  developed  in  this  study  represents  an  entirely  new  effort.  All 
the  properties  mentioned  in  the  physical  description  of  the  upper  atmosphere 
above  have  been  developed  and  incorporated  into  a  computer  program.  The 
sophistication  and  accuracy  of  the  physical  model  are  commensurate  with 
the  uncertainties  in  the  knowledge  of  the  physical  properties  of  the  upper 
atmosphere. 
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The  balance  of  the  report  summarizes  the  research  effort  performed  under 
contract  F19628-69-0258  and  presents  preliminary  results  of  background 
radiance  calculations  for  both  the  lower-  and  upper-atmospheric  regions. 
Section  XI,  for  example,  discusses  the  derivation  of  temperature,  pressure  and 
absorber  profiles  for  the  lower  atmosphere.  These  represent  the  spatial  and 
temporal  variations  of  atmospheric  properties,  and  the  model  atmospheres 
are  used  as  input  data  to  Anding's  radiance  synthesizer  program.  Section 
III  presents  the  results  of  background  radiances  from  the  lower  atmosphere 
in  the  5  to  25  /urn  region,  the  variability,  means  and  extremes  of  infrared  limb 
radiances. 

The  effort  directed  towards  the  upper  atmospheric  portion  of  the  problem  is 
presented  in  Sections  IV  through  VII.  A  discussion  of  infrared  active  gases, 
their  abundances  and  wavelengths  of  emission  is  presented  in  Section  IV. 

These  data  are  input  values  to  the  computational  model  which  is  developed  in 
Section  V.  Preliminary  results  of  infrared  background  radiances  are 
presented  in  Section  VI,  based  upon  the  computational  model.  These  results 
include  examples  of  band  radiances  as  viewed  by  exoatmospheric  and 
endoatmospheric  reference  platforms.  A  detailed  spectra  from  the  I^O 
rotational  band  is  also  included  in  the  results.  The  discussion  in  Section  VII 
concerns  infrared  radiation  from  particulate  material  such  as  clouds  and  dust 
in  the  upper  atmosphere  and  from  dust  in  interplanetary  space.  Consideration 
of  both  gaseous  and  particulate  sources  of  infrared  radiation  provides  a 
reasonably  complete  picture  of  the  natural  infrared  radiant  background  of  the 
near  earth  environment. 
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SECTION  II 

DERIVATION  OF  TEMPERATURE,  PRESSURE,  AND  ABSORBER 
PROFILES  FOR  TIIE  LOWER  ATMOSPHERE  (0  to  70  km) 


Part  of  the  research  effort  during  the  study  was  the  derivation  of  models 
of  the  lower  atmosphere  in  the  altitude  range  from  the  surface  to  70  km. 
These  models  include  the  altitude  distribution  of  temperature,  pressure. 


and  important  absorber  gases  in  the  a  to  25  pi  spectral  region.  The  speci¬ 
fication  of  these  parameters  and  their  distributions  with  altitude  are  used  as 
input  data  to  a  radiation  synthesizer  program  (Anding,  1969)  which  computes 
spectral  radiances  as  a  function  of  arbitrary  view  directions.  The  model 
atmospheres  derived  herein  represent  the  spatial  and  temporal  variations  of 


atmospheric  parameters  about  the  earth.  Extreme  atmospheric  conditions 
are  also  included  as  part  of  the  data  a  t.  sp-  ctral  radiances  computed  from 
these  models  will  give  a  good  indication  of  the  magnitude  and  variation  of  the 
earth's  natural  infrared  background  for  altitudes  up  to  about  70  km.  These 
data  are  important  in  defining  radiometer  design  specifications  and  spatial 
sampling  requirements. 


This  section  discusses  the  sources  of  information  used  to  construct  the  model 
atmospheres,  the  uncertainties  in  our  knowledge  of  parameters,  varia¬ 
tions  as  a  function  of  time,  latitude  and  longitude,  and  the  methods  used  to 
interpolate  parameters  as  a  function  of  altitude.  The  parameter  profiles 
derived  in  the  study  are  tabulated  in  Appendix  A,  and  Appendix  B  summarizes 
the  identification  of  codes  in  A.  Appendix  C  presents  a  selection  of  4  7  model 
atmospheres  that  are  representative  of  spatial  and  temporal  variations,  moan 
and  extreme  conditions  about  the  earth. 


A.  APPROACH  TO  DERIVATION  OF  PROFILES 


The  infrared  radiance  is  related  to  the  atmospheric  temperature,  pressure, 
and  absorber  distributions.  J’o  obtain  estimates  of  the  mean  radiance  and  its 
variability  in  time  and  space  from  theoretical  calculations,  one  must  have 
estimates  of  the  mean  vertical  profiles  of  temperature,  pressure,  and 
absorbers,  and  the  variations  from  the  mean.  The  major  meteorological 
variations  in  the  0  to  70  km  region  are  the  latitudinal  and  seasonal  variations, 
and  season  and  latitude  are  used  as  a  basis  for  deriving  a  set  of  profiles 
representing  atmospheric  variability. 

The  computer  program  used  to  calculate  radiances  for  the  lower  atmosphere 
is  a  modified  version  of  the  one  developed  by  Anding  (1968).  The  absorbing 
gases  included  in  this  program  are  H20  vapor,  CC>2,  Oy  HNOg,  N20  and  CH4. 
Only  the  first  four  of  these  gases  are  incorporated  in  the  5  to  2  5  micron  region. 
Even  though  radiances  from  nitrous  oxide  are  not  computed  in  the  spectral 
region  of  interest  for  the  lower  atmosphere,  profiles  of  N20  are  derived  below 
for  the  sake  of  completeness  and  use  during  future  research  efforts.  Absorption 
constants  for  CH4  are  not  presently  available  for  the  7.  6-micron  band;  however, 
its  omission  should  not  be  significant  in  the  lower  atmosphere  because  of  the 
heavy  absorption  due  to  the  strong  wing  of  the  6.  3-band  of  water  vapor. 


Of  the  five  gases  mentioned,  carbon  dioxide  and  nitrous  oxide  are  uniformly 
mixed  up  to  the  levels  at  which  dissociation  begins.  Water  vapor  and  ozone, 
on  the  other  hand,  vary  vertically,  horizontally,  and  with  time.  Thus,  with 
respect  to  carbon  dioxide  and  nitrous  oxide,  it  is  sufficient  to  specify  a  con¬ 
stant  mixing  ratio  up  to  the  levels  of  dissociation  (for  C02,  ~100  km;  for 
N2(),  ~30  km).  For  both  water  vapor  and  ozone,  however,  as  with  temperature 
and  pressure,  a  number  of  different  profiles  are  required  to  cover  the  range 

of  observed  variations.  Nitric  acid  is  also  variable  with  altitude  and  is 
discussed  in  Appendix  D. 

All  profiles  cover  the  range  from  the  surface  to  90  km.  For  ease  of  tabulation 
and  punching  the  data  on  cards,  and  to  minimize  sources  of  error,  each  profile 
is  initially  tabulated  at  the  minimum  number  of  heights  necessary  to  portray 
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the  main  features  of  the  particular  profile.  These  heights  may  be  called  key 
heights.  Appropriate  interpolation  schemes  are  contained  in  the  computer 
program  that  prepares  the  initial  profiles  for  input  to  the  radiance  calculations. 
The  following  sections  discuss  each  atmospheric  parameter  individually  and 
describe  how  vertical  profiles  representative  of  that  parameter's  mean  state 
and  variability  are  constructed. 


1.  Temperature 


The  U.S.  Standard  Atmosphere  (1962)  temperature  profile  is  adopted  as  the 
mean  temperature  profile.  The  basic  variations  about  the  mean  can  be  rep¬ 
resented  by  several  temperature  profiles  from  the  U.  S.  Standard  Atmosphere 
Supplements  (1966):  The  15°  latitude  annual  mean;  45°  latitude  January  and 
July  profiles;  and  75°  latitude  January  and  July  profiles.  As  presented  in  the 
U.S.  Standard  Atmosphere  Supplements,  the  75°  latitude  profiles  extend  only 
to  30  km.  We  have  extended  these  profiles  to  90  km  as  described  below. 

Observations  of  the  summer  temperature  profile  between  30  and  90  km  at 
Barrow,  Alaska  (71°N)  (Theon,  et  al. ,  1967;  Theon,  1968)  indicate  little 
variability  with  time.  Seven  temperature  profiles  presented  by  Theon  (1968; 
Figure  6)  are  averaged  to  extend  our  75°  July  temperature  profile  to  90  km. 
Although  the  data  are  sparse  for  these  altitudes  in  the  arctic  regions,  the 
available  information  does  indicate  little  variability  in  summer.  Thus,  the 
adopted  profile,  though  based  on  few  observations,  should  be  representative 
of  arctic  summer  conditions  at  these  altitudes. 

For  arctic  (75°  latitude)  winter  conditions,  the  U.S.  Standard  Atmosphere 
Supplements  contains  three  temperature  profiles  from  the  surface  to  3  0  km. 
These  three  —  a  warm,  cold,  and  mean  profile  —  cover  the  range  of  thermal 
regimes  that  exist  in  the  winter  stratosphere.  The  cold  regime  is  dominant 
over  northeastern  Canada  and  the  warm  regime  is  dominant  over  the  Aleutian 
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area.  The  mean  profile  is  extended  to  52.  5  km  using  the  analyses  of  Kantor 
(1966).  Kantor  derives  temperature  at  7  5°N  latitude  between  2  5  and  55  km 
from  radiosonde  data  for  latitude  60^,  meteorological  rocket  wind  and 
temperature  observations  at  Fort  Churchill,  Canada  (55°N)  and  Fort.  Greel.y, 
Alaska  (64°N),  and  application  of  the  thermal  wind  relationship,  l'he  warm 
and  cold  January,  75°  profiles  are  extrapolated  to  52.  5  km  by  assuming  that 
their  temperatures  at  52.  5  km  are  ±  1CTK  from  the  mean  temperature  (see 
Figure  1).  All  three  profiles  are  extrapolated  to  a  common  temperature  of 
210°K  at  90  km,  this  temperature  being  based  upon  the  few  high  latitude  -  high 
altitude  rocket  grenade  observations  (Theon,  et  al,  1967;  Theon,  1968).  It 
should  be  pointed  out  that,  as  opposed  to  summer  conditions  in  the  50  to  90  km 
at  high  latitudes,  the  few  available  soundings  indicate  large  variability  in 
winter  —  up  to  80°K  temperature  change  in  3  hours. 

V 

The  standard  atmosphere,  and  January  temperature  profiles  for  45  N  and  75° N 
(mean)  are  plotted  in  Figure  2;  the  15°  annual  mean,  and  July  temperature 
profiles  for  45°  and  75°  are  plotted  in  Figur  e  3. 

The  temperature  profiles  assembled  thus  far  should  provide  a  reasonable 
estimate  of  the  seasonal  and  latitudinal  variability  of  temperature  between 
the  surface  and  90  km.  In  turn,  they  should  yield  reasonable  estimates  of 
radiance  variability  when  used  as  inputs  to  a  radiance  synthesizer  program. 

Aside  from  the  expected  seasonal  and  latitudinal  variability  of  radiance,  one 
would  also  like  to  estimate  extreme  values  that  might  occur.  For  this  purpose, 
two  extreme  temperature  profiles  have  been  constructed.  These  are  based 
upon  information  on  extreme  temperatures  between  the  surface  and  100,  000 
feet  that  is  presented  in  Valley  (1965,  Figure  2-7).  Hot  temperatures  (those 
exceeded  only  10  percent  of  time  in  the  hottest  regions)  and  cold  temperatures 
(those  exceeded  90  percent  of  the  time  in  the  coldest  regions)  are  presented 
as  a  function  of  altitude.  From  the  average  range  between  hot  and  cold 
temperatures  over  this  altitude  span  which,  for  a  normal  distribution  represents 
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Figure  1.  Adopted  temperature  profiles  for  January,  75°  Latitude 
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figure  2.  Standard  Atmosphere,  45°  January,  and  75°  January 
(mean)  Temperature  Profiles 
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2.6  a  standard  deviation  of  16  K  is  derived.  For  the  extreme  temperature 
profiles,  a  3  a  level,  or  plus  or  minus  ~5CPK  is  adopted.  By  adding  50°K  to  the 
standard  atmosphere  temperature  profile  at  all  levels,  one  obtains  the  extreme 
hot  profile;  by  subtracting  50°K,  one  obtains  the  extreme  cold  temperature 
profile.  These  synthetic  temperature  profiles  should  be  useful  for  estimating 
extremes  in  infrared  radiance. 

Temperatures  at  key  heights  from  the  surface  to  90  km  for  all  the  temperature 
profiles  are  tabulated  in  Appendix  A.  Linear  interpolation  may  be  used 
between  key  heights. 


2.  Pressure 

Pressure  profiles  are  calculated  from  the  temperature  profiles  with  the  use 
of  the  hydrostatic  equation.  A  surface  pressure  appropriate  to  each  of  the 
Supplemental  Atmospheres  is  used.  For  the  two  extreme  temperature  pro¬ 
files,  a  surface  pressure  of  1013  mbar  is  used.  The  following  integrated  forms 
of  the  hydrostatic  equation  are  used  in  the  pressure  computations: 

For  constant  lapse-rate  layers: 


P(z)  =  P(z  ) 
o 


T(z) 


g/ctR 


T<VJ 


For  isothermal  layers: 

P(z)  =  P(z  )  exp 


g(z  -  zq) 

R  T(z  ) 
o 


(1) 


(2) 


where  P  is  pressure,  z  is  altitude,  z^  is  the  altitude  of  the  base  of  a  layer, 
T  is  temperature,  g  is  the  gravitational  acceleration,  a  is  the  temperature 
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lapse  rate  in  the  layer,  and  R  is  the  gas  constant  for  air.  The  variation  of 
the  gravitational  acceleration  g  with  altitude  is  included  in  the  pressure 
computation. 


3.  Water  Vapor 

For  purposes  of  constructing  profiles  representative  of  the  variability  of 
water  vapor,  the  0  to  90  km  altitude  range  is  divided  into  two  regions:  a  lower 
atmospheric  region  extending  from  the  surface  to  1  5  km  (troposphere),  and 
an  upper  atmospheric  region  extending  from  15  to  90  km  (which,  in  this 
section,  is  referred  to  as  the  stratosphere).  The  reason  for  this  division 
is  the  fact  that  there  is  very  little  information  on  the  water  vapor  content  of 
the  stratosphere  and  quite  a  large  amount  of  data  on  the  water  vapor  content 
of  the  troposphere  -  especially  the  lower  troposphere.  Thus,  IJ00  profiles 
in  the  troposphere  are  representative  of  seasonal  and  latitudinal  vuriuijil 
for  the  stratosphere,  only  estimates  of  the  mean  HgO  profile  and  p.  ^sibJc 
extremes  are  possible. 

Mean  vertical  profiles  of  atmospheric  moisture  content  have  been  prepared 
(Valley,  1965)  from  the  surface  to  10  km  from  the  Supplemental  Atmospheres. 
These  data  are  used  as  the  basis  for  the  estimates  of  the  variability  n  iihiu  a 
season  and  latitude  belt,  the  information  provided  in  the  Atmospheric  Humidity 
Atlas  —  Northern  Hemisphere  (Gringorten,  et  al. ,  1966)  is  used.  This 
publication  contains  hemispheric  distribution  of  various  percentile  values  of 
water  vapor  content  for  the  surface  and  for  the  850-,  700-,  500-,  and  400-mbar 
levels  as  a  function  of  month  of  the  year.  The  analyses  are  based  upon  five 
years  of  world  wide  radiosonde  data.  From  these  maps,  estimates  of  lr'rh 
(95  percentile)  and  low  (5  percentile)  water  vapor  mixing  ratios  for  seasons 
and  latitudes  are  obtained,  similar  to  those  used  for  the  temperature  profiles 
and  mean  water  vapor  profiles. 


When  the  water  vapor  mixing  ratio  is  plotted  on  a  semilogarithmic  graph  versus 
alMtude,  it  can  be  seen  that  the  logarithm  of  the  mixing  ratio  decreases  ap¬ 
proximately  linearly  with  altitude.  Using  the  data  derived  from  the  sources 
listed  above,  the  high,  mean  and  low  H90  profiles  are  derived  as  illustrated 
in  1H  igure  4.  The  plotted  points  are  for  45  N,  July,  and  are  from  the  original 
sources;  the  solid  lines  represent  the  adopted  distributions;  and  the  short 
horizontal  bars  divide  the  troposphere  into  layers  in  which  the  logarithm  of 
the  water  vapor  mixing  ratio  decreases  at  a  constant  linear  rate  with  altitude. 
Profiles  were  constructed  in  a  similar  manner  for  15°  annual,  45°  .January, 
and  75°  January  and  July  conditions. 

1  or  the  oveiall  mean,  or  standard  atmosphere,  water  vapor  profile  in  the 
troposphere,  the  mixing  ratio  profile  presented  by  Sissenwine,  et  al.  (1968) 
is  adopted  which  is  based  upon  the  data  in  Gringorten,  et  al.  (1966). 

Because  of  the  low  mixing  ratios  and  pressures,  water  vapor  in  the  strato¬ 
sphere  is  extremely  difficult  to  measure.  Standard  radiosonde  sensors 
generally  cannot  measure  water  vapor  concentrations  at  pressures  much  below 
400  mbar.  A  recent  summary  of  available  measurements  has  been  presented 
in  Leovy  (1969)  and  is  quoted  below. 

For  observations  in  which  serious  contamination  effects  seem 
to  have  been  avoided,  both  hygrometer  and  spectroscopic  obser¬ 
vations  indicate  remarkably  uniform  water  vapor  mixing  ratios 
in  the  range  1-4  x  10  6  gm/gm  from  just  above  the  tropopause 
up  to  20  or  2  5  km.  Above  2  5  km,  the  picture  is  not  so  clear. 

Relatively  few  measurements  have  been  made  at  these  heights. 

Spectroscopic  measurements  are  subject  to  ambiguity  of  inter¬ 
pretation  above  20  km  (Houghton,  1960;  and  Williamson,  1965). 

Some  frost-point  measurements  have  indicated  continued  dryness 
with  mixing  ratios  of  1-4  x  10  6  gm/gm  up  to  30  km  (Goldsmith, 

1J64).  In  a  year-long  series  of  frost  point  hygrometer  ascents  at 
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rofiles  Indicated  by  Lines. 


Washington,  D.  C.,  and  Trinidad,  Mastenbrook  (1968)  has  found 
mixing  ratios  predominantly  in  the  range  1-4  x  1  0"  gm/gm  above 
20  km,  but  with  some  seasonal  variations,  and  with  greater  moisture 
concentrations  at  the  higher  levels  on  a  few  ascents.  Brown,  et  al. 

(1961)  flew  a  cooled  vapor  trap  on  a  balloon  between  25  and  28  km 
and  found  mixing  ratios  ranging  from  9  to  77  x  10  gm/gm  at  these 
altitudes.  They  were  able  to  rule  out  the  possibility  of  contamination 
from  the  balloon  by  introducing  a  small  amount  of  deuterium  into  the 
balloon  as  a  tracer.  The  possibility  remains,  however,  of  some  con¬ 
tamination  from  other  parts  of  their  system,  although  great  care  was 
taken  to  eliminate  all  sources  of  contamination.  Recent  frostpoint 
measurerm  nts  by  Sissenwine,  et  al.  (1968),  also  indicate  increasing 

water  vap>  r  concentration  above  20  km,  with  mixing  ratios  of 

_  0 

20  x  10  gm/gm  at  25  km.  The  occasional  occurrence  of  mother-of- 
pearl  clouds  in  this  portion  of  the  stratosphere  has  been  cited  as 
indirect  evidence  for  mixing  ratios  as  high  as  10~4  gm/gm  (Newell, 

1963)." 

As  the  above  discussion  indicates,  there  is,  at  present,  no  consensus  as  to  the 
mean  H20  mixing  ratio  profile  in  the  stratosphere  nor  of  its  variability.  For 
the  present  purpose  of  providing  inputs  to  radiance  calculations,  the  three 
H20  profiles  for  the  stratosphere  are  adopted.  These  should  cover  the  range 
of  present  uncertainty  in  this  gas  concentration  as  well  as  any  natural 
variability. 

For  the  mean  stratospheric  water  vapor  profile,  the  mid-latitude  moisture 
profile  presented  by  Sissenwine  et  al.  (1968)  are  adopted.  Up  to  32  km  this 
profile  is  based  upon  mid-latitude  stratospheric  frost  point  observations 
obtained  by  Sissenwine  and  his  colleagues.  Between  32  km  and  80  km,  it  is 
based  upon  an  extrapolation  to  a  dewpoint  of  -130°C  at  80  km.  This  dewpoint 
is  based  upon  an  observation  of  the  occurrence  and  disappearance  of  noctilucent 
clouds  (Anonymous,  1963).  The  adopted  profile  has  a  maximum  mixing  ratio 
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of  19  ppm  at  26  km.  The  mixing  ratio  decreases  with  altitude  above  this  level 
to  a  value  of  0.  6  ppm  at  80  km.  Between  80  and  90  km,  a  constant  mixing 
ratio  is  assumed. 

Our  dry  stratospheric  water  vapor  profile  is  based  mainly  upon  the  frost - 
point  observations  of  Mastenbrook  (1968).  A  value  of  1  ppm  between  15  and 
30  km  is  assumed,  a  drop  to  0.  1  ppm  at  80  km,  and  a  constant  value  between 
80  and  90  km. 

For  the  wet  stratosphere,  a  value  of  50  ppm  between  15  and  30  km  is  assumed, 
a  drop  to  5  ppm  at  80  km,  and  a  constant  value  between  80  and  90  km. 

The  adopted  standard  atmosphere  water  vapor  profile  for  the  troposphere  and 
the  dry,  mean,  and  wet  stratospheric  water  vapor  profiles  are  illustrated  in 
Figure  5.  Just  as  the  wet,  mean,  and  dry  stratospheric  profiles  can  be  placed 
on  top  of  the  standard  atmosphere  tropospheric  water  vapor  profile,  so  can 
they  be  connected  to  any  of  the  seasonal  and  latitudinal  tropospheric  water 
vapor  profiles.  From  the  top  of  tropospheric  H^O  profiles  (10  km)  to  the  base 
of  the  stratospheric  H^O  profiles  (15  km)  a  simple  linear  change  with  altitude 
of  the  logarithm  of  the  mixing  ratio  is  assumed. 

Water  vapor  mixing  ratios  at  key  altitudes  between  the  surface  and  10  km  for 
all  the  tropospheric  profiles,  and  between  15  and  90  km  for  the  three  strato¬ 
spheric  profiles,  are  tabulated  in  Appendix  A.  Interpolation  between  key 
heights  is  accomplished  by  assuming  that  the  logarithm  of  the  mixing  ratio  is 
a  linear  function  of  altitude. 


4.  Ozone 

Since  1963,  a  small  network  of  North  American  stations  has  been  flying 
chemical  ozonesondes  to  measure  the  vertical  distribution  of  atmospheric 
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Figure  5 


.  Adopted  Standard  Atmosphere  Water  Vapor  Profile 
for  Troposphere  and  Dry,  Mean,  and  Wet  Stratosphe 
Water  Vapor  Profiles 
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ozone  between  the  surface  and  balloon  burst  altitudes  (~  3  0  km).  An  interim 
summary  of  these  data,  covering  the  year  s  1963  to  1964,  has  been  published 
by  Hering  and  Borden  (1965).  For  constructing  ozone  profiles  to  3  5  km  in 
the  present  study,  we  have  used  unpublished  ozone  cross-sections  recently 
prepared  by  Hering  and  Borden  (1969)  and  covering  all  of  the  ozone  soundings 
to  date.  Between  3  5  and  90  km,  where  ozone  observations  are  scarce,  the 
distributions  suggested  by  Kuhn  (1966)  are  used. 

The  Hering-Borden  cross-sections  extend  from  equator  to  north  pole  and 
from  the  surface  to  10  mbar.  For  each  season,  there  is  one  cross-section 
for  ozone  density  and  one  for  the  standard  deviation  of  the  ozone  density.  In 
order  to  keep  the  number  of  model  atmospheres  within  a  reasonable  number 
(47)  and  to  be  consistent  with  the  temperature  and  water  vapor  profiles,  the 
ozone  data  are  grouped  together  into  two  seasons,  represented  as  winter  and 
summer.  Based  on  the  ozone  density  cross-sections,  the  vertical  ozone  pro¬ 
files  for  15°  annual  mean,  4  5°  winter  mean,  4  5°  summer  mean,  75°  winter 
mean,  and  75  summer  mean  are  constructed.  Based  on  the  standard  deviation 
cross-sections,  the  extreme  profiles  for  each  of  the  latitude-season  zones 
listed  above  are  derived.  These  extreme  profiles  are  constructed  by  adding 
or  subtracting  a  2cr  variability  to  or  from  the  mean  profiles  for  each  of  the 
latitutde-season  zones.  Both  the  mean  profile  and  the  extreme  ozone  were 
extended  graphically  from  ~30  to  35  km.  The  spring  and  fall  cross-sections 
given  by  Hering-Borden  are  similar  and  fall  within  the  models  indicated  above. 


For  our  overall  mean,  or  standard  atmosphere,  ozone  profile  in  the  lowest 
3  5  km,  the  standard  atmosphere  ozone  profile  presented  in  the  U.  S.  Standard 
Atmosphere  Supplements  (1966)  is  adopted,  which  is  based  upon  the  results  of 
ozonesonde  flights  at  mid-latitude  stations. 

Between  35  and  90  km  the  shape  of  the  ozone  profile  suggested  by  Kuhn  (1966) 
is  used.  Kuhn  assumes  that  the  ozone  mixing  ratio  is  proportional  to  the 
pressure  raised  to  a  power.  This  may  be  written  as: 
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03  <z>  =  °3  (Zj)  j^P(z)/P(Z1)J  (3) 

where  03  is  the  ozone  mixing  ratio,  z  is  altitude,  Zj  is  a  reference  altitude, 
and  P  is  pressure.  Between  3  5  and  50  km,  p  is  based  mainly  on  photochemical 
equilibrium  calculations  and  is  equal  to  0.46.  Between  50  and  90  km,  p  is 
based  upon  rocket  observations,  which  agree  with  photochemical  calculations, 
and  is  equal  to  0.  614.  Photochemical  theory  suggests  large  dirunal  variations 
in  the  upper  layers  of  the  ozone  profile  (>  50  km).  These  are  not  presently 
included  in  the  ozone  profiles,  but  can  easily  be  programmed  since  the  region 
between  3  5  and  90  km  can  be  subdivided  into  a  number  of  layers  with  different 
reference  levels  z,  and  values  of  p.  Thus,  as  opposed  to  the  0  to  3  5  km 
region,  where  the  C>3  mixing  ratios  are  specified  at  key  heights  and  the  O 
profile  for  any  desired  vertical  resolution  is  obtained  by  linear  interpolation, 
the  Os  profile  above  3  5  km  is  calculated  analytically  by  means  of  a  group  of 
formulas  similar  to  equation  (3). 

Figure  6  illustrates  the  adopted  standard  atmosphere  ozone  profile.  The  ozone 
mixing  ratios  at  key  heights  between  the  surface  and  3  5  km  for  all  the  ozone 
profiles  are  tabulated  in  Appendix  A.  Mixing  ratios  between  key  heights  are 
obtained  by  linear  interpolation  in  altitude. 


5.  Nitrous  Oxide 

Bates  and  Hays  (1967)  have  recently  reviewed  information  on  nitrous  oxide 
and  performed  theoretical  calculations  of  vertical  profiles.  The  most  probable 
source  of  atmospheric  N^O  is  microbiological  action  at  the  earth's  surface. 

The  N2()  produced  at  the  earth's  surface  is  mixed  by  turbulent  diffusion  through¬ 
out  the  troposphere  and  into  the  stratosphere.  A  mean  tropospheric  fractional 
volume  abundance  of  2.  5  x  10"7  is  suggested  by  the  few  available  observations. 

In  the  stratosphere,  photodissociation  takes  place  and  the  N20  abundance 
probably  decreases  with  altitude. 
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Figure  6.  Adopted  Standard  Atmosphere  Ozone  Prof! 


Bates  and  Hays  (1967)  have  computed  vertical  profiles  of  f(N90),  the  fractional 
\  olume  abundance,  from  the  surface  to  40  km,  based  upon  a  simple  model  in 
which  vertical  eddy  diffusion,  a  surface  source,  and  an  atmospheric  photo¬ 
dissociation  sink  are  considered.  The  resulting  profiles  which  are  plotted  in 
Figure  7  indicate  that  the  stratospheric  profiles  are  sensitive  to  the  value  of 
*  the  stratospheric  eddy  diffusion  coefficient. 


These  computed  N20  profiles  have  been  extended  to  90  km  by  noting  that  log 
f(N2( B  in  the  range  1  or  2  x  10  to  10  decreases  with  altitude  at  roughly 
the  same  rate,  independent  of  K^.  Thus,  for  heights  above  the  level  at  which 
f(N2Q)  =  10  ,  the  following  formula  is  used  to  represent  the  N  O  concentration: 


N2o 


(z)  = 


fN00(V 


■p(z  - 


V 


(4) 


where 

p  -  0.  24  km  * 

the  adopted  N20  profiles  are  plotted  in  Figure  8  for  four  possible  values  of 

the  mean  stratospheric  eddy  diffusion  coefficient.  The  adopted  N20  profile 

lor  the  standard  atmosphere  incorporates  a  diffusion  coefficient  of  K  = 

4  2-1  z 

10  cm  sec  .  The  adopted  N20  profiles,  from  the  surface  to  40  km,  are 

tabulated  in  Appendix  A.  Above  40  km,  equation  (4)  is  used. 

It  should  be  noted  again  that  the  computer  program  which  calculates  radiances 
for  the  lower  atmosphere  does  not  include  nitrous  oxide  in  the  5  to  2  5  /um 
spectral  region.  However,  the  program  does  include  this  gas  at  wavelengths 
less  than  5  gm.  Therefore,  N2Q  profiles  were  derived  during  this  study  for 
the  sake  of  completeness  and  use  during  future  research  efforts. 
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Figure  7.  Fractional  Volume  Abundance  of  Atmospheric 
Nitrous  Oxide,  f(N20),  for  Selected  Values  of 
the  Vertical  Exchange  Coefficient  Kz  (after 
Bates  and  Hays,  1967) 
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Figure 


6.  Carbon  Dioxide 


For  our  present  purposes,  carbon  dioxide  may  be  assumed  to  be  uniformly 
mixed  from  the  surface  to  ~  100  km.  Recent  aircraft  observations  (Georgii 
and  Jost,  1969)  suggest  a  decrease  from  the  troposphere  to  the  stratosphere 
of  the  order  of  1  percent  of  mean  mixing  ratio.  However,  such  a  small 
variation  is  of  little  concern  to  the  current  problem.  There  is  some  un¬ 
certainty  in  the  value  of  the  mean  atmospheric  CC>2  mixing  ratio.  Since  the 
beginning  of  the  century,  there  has  been  a  10  to  15  percent  increase  in 
atmospheric  CC>2  (Callendar,  1958)  as  a  result  of  increased  burning  of  fossil 
fuels.  Recent  studies  (Bolin  and  Keeling,  1960;  Keeling,  1960)  suggest 
that  this  increase  is  continuing  at  a  rate  of  the  order  of  1  ppm  per  year.  A 
recent  survey  of  atmospheric  carbon  dioxide  variability  (Peterson  et  al. , 

1967)  suggests  a  mean  C02  mixing  ratio  of  314  ppm  by  volume  applicable  some 

four  to  six  years  ago.  A  CC>2  mixing  ratio  of  320  ppm  by  volume  is  adopted 
in  this  study. 


7.  Nitric  Acid 


Nitric  acid  has  been  identified  as  a  minor  constituent  in  the  atmosphere 
(Murcray,  et  al,  1968)  with  a  strong  absorption  band  near  11.2  microns.  The 
computer  program  developed  by  Anding  (1968)  does  not  include  a  model  for 
nitric  acid.  Calculated  radiances  in  the  window  region  from  10  to  12  microns 
have  shown  extremely  low  values;  therefore,  it  has  been  necessary  to 
synthesize  a  model  for  absorption  by  nitric  acid.  The  inclusion  of  nitric  acid 
represents  the  major  change  to  the  University  of  Michigan  program. 


Transmission  by  nitric^acid  is  given  by  the  following  empirical  relationship: 


t  =  exp 


0. 84  x  10 


-12 


(4a) 
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atmospheric  pressure 

pressure  at  normal  conditions  (760  mm  Hg) 
atmospheric  temperature 

temperature  at  normal  conditions  (288.16°K) 
altitude 

path  distance  along  line  of  sight 


The  development  of  the  model  is  given  in  Appendix  D. 


SECTION  III 

VARIABILITY  OF  INFRARED  LIMB  RADIANCES 
IN  THE  LOWER  ATMOSPHERE 


Infrared  limb  radiances  were  computed  for  a  selected  sample  of  model  at¬ 
mospheres  derived  in  Section  II.  These  radiances  represent  the  mean  con¬ 
ditions,  the  extreme  conditions  and  the  uncertainties  in  radiance  magnitudes 
due  to  our  knowledge  of  atmospheric  parameters  in  the  stratosphere  and  the 
mesophere.  The  results  presented  herein  include  both  spectral  radiances  and 
band  radiances  for  absorbing  gases,  computed  as  a  function  of  tangent  height. 
Both  spatial  and  temporal  variations  of  band  radiances  are  discussed,  and 
extreme  conditions  are  presented  as  a  function  of  tangent  height.  These  re¬ 
sults  are  useful  in  specifying  radiometer  requirements  and  defining  data 
sampling  requirements  for  tangent  heights  up  to  about  70  km. 

For  all  calculations,  carbon  dioxide  was  assumed  to  be  uniformly  mixed 
throughout  the  atmosphere  with  a  volumetric  mixing  ratio  of  320  ppm.  The 
particular  meteorological  profile  used  in  a  radiance  calculation  is  represented 
by  a  four  number  code,  the  first  number  indicating  the  type  of  temperature 
profile,  the  second  number  indicating  the  type  of  trophospheric  water  vapor 
profile,  the  third  number  indicating  the  type  of  stratospheric  water  vapor 
profile,  and  the  fourth  number  indicating  the  type  of  ozone  profile.  The 
codes  are  explained  in  Appendix  B.  For  example,  atmospheric  profile 
10-12-3-20  represents  a  hot  extreme  temperature  profile,  a  tropospheric 
water  vapor  profile  typical  of  an  extremely  moist  condition  at  15°  latitude 
(this  is  the  wettest  tropospheric  water  vapor  profile  in  our  data  set),  a  wet 
stratospheric  water  vapor  profile,  and  an  ozone  profile  typical  of  an  extremely 
high  ozone  content  at  75°  latitude  during  the  winter  season  (this  is  the 
ozone  profile  with  the  highest  total  amount  in  our  data  set). 
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F  *6ure  9  illustrates  the  5  to  25  Mm  spectrum  of  the  standard  atmosphere 
(code  1-1-1 -1)  for  three  different  tangent  heights:  h  =  50  km,  60  km,  and 
70  km.  The  major  features  of  the  spectra  are  the  6.3  Mm  H20  band,  the 
9.6  Mm  Os  band,  the  15  Mm  CC2  band,  and  the  rotational  H20  band  beyond 
18  Mm.  Ihese  features  are  seen  in  emission  since  they  are  being  viewed 
against  the  cold  background  of  space. 

The  general  decrease  of  radiance  with  increasing  tangent  height  in  the  50  to 
70  km  range  is  due  to  decreasing  absorber  amounts  and  temperatures  with 
altitude.  Highest  emissions  are  in  the  middle  of  the  15  Mm  CC>2  band;  lowest 
emissions  are  in  the  window  region  near  10  Mm.  It  is  noted  that  the  range 
in  radiance  between  these  two  regions  covers  almost  five  orders  of  magnitude. 

Figure  10  shows  how  the  spectrum  for  tangent  height  60  km  might  vary  under 
extreme  atmospheric  conditions.  The  topmost  spectrum  represents  extremely 
hot,  wet  and  high  ozone  conditions  (10-12-3-20);  the  middle  spectrum  repre¬ 
sents  the  standard  atmosphere  (l-l-l-l);  and  the  lowest  spectrum  represents 
extremely  cold,  dry,  and  low  ozone  conditions  (9-19-2-17). 

The  range  between  the  top  and  bottom  spectra  is  probably  representative  of 
the  maximum  range  to  be  expected  at  this  tangent  height,  and  the  actual 
magnitudes  are  estimates  of  the  highest  and  lowest  spectral  radiances  to 
be  expected  at  this  tangent  height. 

One  of  the  major  uncertainties  in  the  calculations  discussed  above  is  the  water 
vapor  distribution  in  the  stratosphere.  To  illustrate  the  effect  of  this  uncer¬ 
tainty  in  the  spectra  in  the  6.3m  HgO  band  and  rotational  HsO  band,  computa¬ 
tions  were  performed  with  a  standard  atmosphere  temperature  profile  with  a 
dry  and  wet  stratosphere.  The  computed  spectra  for  the  two  water  vapor 
bands  for  a  tangent  height  of  60  km  are  plotted  in  Figure  11  along  with  the 
spectra  for  the  adopted  mean  stratospheric  H20  profile.  It  can  be  seen  that 
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Figure  9.  Infrared  Spectra  at  Tangent  Heights  of  50,  60  and  70  km 
for  Standard  Atmosphere  Temperature,  Water  Vapor, 
and  Ozone  Conditions 
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Figure  10.  Infrared  Spectra  at  a  Tangent  Height  of  60  km  for 
Standard  Atmosphere  Conditions  and  Extreme 
Conditions  on  Temperature,  Water  Vapor,  and  Ozone 
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Figure  11.  Infrared  Spectra  at  a  Tangent  Height  of  60  km 
for  the  H2O  Bands  for  Different  Stratospheric 
H2O  Mixing  Ratios 
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the  range  in  computed  spectral  radiance  between  dry  (1-19-2-20)  and  wet 
(1-12-3-20)  stratospheres  can  be  almost  a  full  order  of  magnitude. 

Of  particular  interest  is  the  variation  of  radiance  with  tangent  height  (horizon 
radiance  profile)  in  the  absorption  bands  and  in  the  windows.  This  is  illus¬ 
trated  in  Figures  12  through  14  for  the  6.  3  p m  HgO  band  (5.  3  pm  to  7.  3  gm), 
the  9.  6  gm  band  (9.  3  gm  to  10.  3  gm),  and  the  rotational  band  of  HgO 
(18.0  gm  to  25  gm).  Each  figure  contains  three  horizon  radiance  profiles. 
Each  horizon  radiance  profile  is  computed  with  the  standard  atmosphere 
temperature  profile  but  with  different  absorber  (HgO  and  O^)  amounts  ranging 
from  high  (1-12-3-20)  to  mean  (1-1-1  —  1)  to  low  (1-19-2-17)  absorber  amounts. 
Thus,  on  one  graph  the  differences  between  horizon  radiance  profiles  are  due 
solely  to  differences  in  absorber  amount  for  that  particular  absorption  band. 

The  general  decrease  of  radiance  with  tangent  height  is  due  mainly  to  the 
general  decrease  of  absorber  concentration  with  altitude,  and  also,  for  seme 
parts  of  the  atmosphere,  to  a  decrease  of  temperature  with  altitude.  In  the 
9.6  gm  Og  band,  there  is  actually  a  slight  increase  of  radiance  with  tangent 
height  in  the  vicinity  of  20  km  tangent  height.  This  is  due  to  the  increase  of 
both  ozone  concentration  and  temperature  in  this  region  of  the  atmosphere. 

In  most  cases,  the  highest  radiance  at  a  particular  tangent  height  is  associated 
with  the  highest  absorber  concentration.  However,  at  the  low  tangent  height  of 
the  horizon  radiance  profiles  for  the  water  vapor  bands,  this  situation  is  re¬ 
versed.  The  highest  radiances  are  associated  with  the  lowest  water  vapor 
concentrations.  This  reversal  occurs  because  the  model  with  the  highest  HgO 
concentration  reaches  opacity  at  some  tangent  height  above  the  surface.  For 
this  model,  any  further  decrease  of  tangent  height  will  not  result  in  substantial 
increase  in  radiance  because  only  a  small  fraction  of  the  emitted  radiation 
from  these  lower  layers  escapes  to  space.  For  the  model  with  the  lowest 
HgO  concentration,  opacity  is  not  reached  and  the  radiance  increases  with 
decreasing  tangent  height  along  with  the  temperature  up  to  the  zero  tangent 
height  (surface)  level. 
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Band  RADIANCE  (WATTS/CM2  STtR) 


Figure  12.  Horizon  Radiance  Profiles  for  the  6.3  /i m 
H2O  Band  for  Standard  and  Extreme  HoO 
Profiles 


33 


40  50  60  70 

IEIGHT,  h  (Km) 


lies  for  the  9.  6  jum  O 
Extreme  O3  Profiles 


TANGENT  HEIGHT,  h  (Km) 


Figure  14.  Horizon  Radiance  Profiles  for  the  18-2  5  jum  Region 
for  Standard  and  Extreme  HnO  Profiles 


Examination  of  the  horizon  radiance  profiles  for  the  9.  6  gm  O  band  reveals 
that,  above  a  tangent  height  of  10  km  the  radiances  are  highest  for  the  stan¬ 
dard  atmosphere  Og  profile  rather  than  for  our  adopted  high  C>3  amount  pro- 
flle*  t  he  reason  for  this  is  that  the  high  Og  abundance  profile  is  based  upon 
an  extreme  high  C>3  profile  that  might  occur  at  high  latitudes  in  winter 
l  he  shaPe  of  this  profile  is  such  that  much  of  the  ozone  concentrated  in  the 

iavt  I  below  ’  1  klT  *  rhis  sliaPe  is  different  from  the  standard  atmosphere 
profile,  which  has  a  higher  center  of  gravity.  As  a  result,  Og  mixing  ratios 
above  an  altitude  of  about  28  km  are  actually  higher  in  the  standard  atmosphere 
profile  rather  than  in  the  high  Og  profile,  although  the  total  amount  of  O  is 
hiaht  st  in  the  high  Og  model.  This  difference  in  shape  combined  with  the 
absorption  characteristics  of  the  9.6  gm  band  causes  the  highest  radiances 

above  10  km  tangent  height  to  be  associated  with  the  standard  atmosphere  O. 
distribution.  3 


The  rapid  fall-off  of  radiance  with  increasing  tangent  height  characteristic  of 
the  absorption  bands  at  the  high  tangent  heights  is  found  throughout  the  entire 
horizon  radiance  profile  in  the  window  region  at  10.  3  gm  (Figure  15).  Over 
the  range  of  tangent  heights  zero  to  60  km,  the  10.  3  gm  radiance  drops  almost 
seven  orders  of  magnitude.  For  the  window  computations,  besides  varying 
the  water  vapor  profile  through  its  extremes,  the  temperature  profile  was 
simultaneously  varied  through  its  extremes.  Thus,  the  difference  between 
top  and  bottom  horizon  radiance  profiles  in  Figure  15  represents  the  extreme 
range  over  which  the  10.3  gm  radiance  might  vary. 


Figures  16  through  19  illustrate  the  extreme  ranges  of  radiance  to  be  ex- 

ThC  L  '!  *he  abS°rPti0n  banC'S  a*  tangent  hei§hts  °f  ze™,  20,  40,  and  60  km. 

.  "  ladia-‘ceb  represent  integrals  over  the  bandwidth  indicated  by  the 

uar  mes.  The  range  in  band  radiance  of  the  6.3  gm  H  O  band  is  consistently 
gi  eater  than  the  range  of  rotational  band  due  to  the  greater  sensitivity  of 
la  nek]  an  emission  to  temperature,  at  the  lower  wavelengths.  These  are 
ex  reme  ranges  and  should  be  useful  in  planning  measurement  programs. 
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Figure  15,  Horizon  Radiance  Profiles  for  the  10,3  p m  window  for 
Standard  and  Extreme  HgO  and  Temperature  Profiles 
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WAVELENGTH  X  (pM) 


Figure  16.  Extreme  Ranges  of  Band  Radianc 
for  Tangent  Height  0  km 


Figure  17.  Extreme  Ranges  of  Band  Radiance  for  Tangent  Height  20  km 
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WAVELENGTH  X  (MM) 


Figure  18.  Extreme  Ranges  of  Band  Radiance  for  Tangent  Height  40  km 


•  • 


WAVELENGTH  X  (hM) 


Figure  19.  Extreme  Ranges  of  Band  Radiance  for  Tangent  Height  60  km 


Also  useful  in  planning  measurement  programs  is  the  average  spatial  varia¬ 
tion  of  radiance  over  the  globe,  figures  20  and  21  show  the  mean  latitudinal 
variations  of  band  radiance  at  30  km  tangent  height  for  winter  and  summer, 
respectively.  The  decrease  of  radiance  toward  the  pole  during  winter  and  the 
slight  increase  toward  the  pole  during  summer  reflects  mainly  the  upper 
stratospheric  temperature  gradient  during  these  seasons. 

Figure  22  shows  the  variability  in  band  radiance  at  30  km  tangent  height  that 
might  be  experienced  as  a  result  of  longitudinal  temperature  gradients  in  the 
arctic  winter  stratosphere.  These  temperature  gradients  arise  as  a  result 
of  the  generally  warm  Aleutian  stratosphere  and  cold  northeastern  c  anada 
stratosphere  during  winter. 

The  sample  of  calculations  discussed  above  was  selected  to  obtain  informa¬ 
tion  on  the  mean  spectra  and  the  variability  about  the  mean.  In  particular, 
calculations  were  performed  with  meteorological  profiles  representing  ex¬ 
tremes  of  temperature  and  absorber  concentrations.  These  should  provide 
valuable  information  on  the  high  and  low  radiances  to  be  expected.  Calcu¬ 
lations  were  also  performed  with  meteorological  profiles  representative  of 
average  spatial  variations.  These  should  provide  information  in  the  normal 
horizontal  variability  of  infrared tiffiti  radiance. 

Spectral  radiances  for  47  atmospheric  models  are  given  in  Appendix  E  in  10- 
kilometer  steps.  The  i.  jdel  numbers  and  tabulated  profile  parameters  are 
given  in  Appendix  C 
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LATITUDE  (DEG) 


Figure  20. 


Mean  Latitudinal  Variation  of  Band  Radiances  at 
h  =  30  km  During  Winter 


SECTION  IV 


INFRARED  ACTIVE  GASES  IN  THE 
UPPER  ATMOSPHERE 

A.  GENERAL  CONSIDERATIONS 

Knowledge  of  the  day  and  night  airglow  emissions  in  the  infrared  region  of 
the  electromagnetic  spectrum  is  of  great  importance  to  military  surveillance 
schemes.  In  addition,  measurement  of  such  radiation  can  provide  the 
scientist  with  information  on  the  energy  content  of  the  atmosphere,  its 
chemical  composition  and  the  physical  reactions  which  take  place. 

All  molecules  except  homonuclear  diatomic  molecules  such  as  N2,  09  and 
II ^  are  potential  strong  emitters  of  infrared  radiation.  This  section  dis¬ 
cusses  estimates  of  abundances  of  infrared  active  molecules  and  describes 
the  considerations  which  led  to  the  selection  of  CC>2,  H20,  03,  NO,  NgO  and 
CH^  alone  for  inclusion  in  preliminary  high  altitude  5  to  25  Mm  infrared 

radiance  calculations.  Since  practically  all  infrared  vibration-rotation 

-2  - 1 

bands  have  strengths  of  the  order  of  100  to  1000  cm  ATM  ,  it  was 
possible  to  make  this  selection  on  the  basis  of  chemical  abundances  and  the 
position  of  infrared  vibration-rotation  bands. 

Most  of  the  infrared  emission  of  the  earth's  atmosphere  comes  from  mole¬ 
cules  present  as  minor  constituents.  The  energy  of  an  infrared  quantum  of 
radiation  is  such  that  in  general  it  corresponds  to  the  difference  in  energies 
of  vibrational  or  rotational  levels  within  a  single  electronic  state  of  a  mole¬ 
cule.  In  order  for  a  transition  resulting  in  emission  of  radiation  to  take 
place  between  two  such  levels,  there  must  be  a  change  in  the  electric  dipole 
moment  of  the  molecule.  This  is  not  possible  in  the  case  of  homonuclear 
diatomic  molecules  such  as  nitrogen  or  oxygen.  The  emission  spectra  of 
these  molecules  arises  from  transitions  between  electronic  states,  and 
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under  normal  conditions  their  spectra  does  not  extend  to  wavelengths  as 
long  as  5  The  electric  dipole  moment  of  heteronuclear  molecules  such 

as  nitric  oxide  (NO)  and  carbon  monoxide  (CO)  is  different  for  different 
vibrational  levels,  and  they  have  prominent  infrared  vibrational  spectra. 

All  polyatomic  molecules  such  as  carbon  dioxide  (COg)  have  infrared  vibra¬ 
tional  spectra.  In  addition,  molecules  with  a  permanent  electric  dipole 
moment  such  as  CO,  NO  and  water  vapor  (HgO)  have  pure  rotational  spectra 
in  the  far  infrared. 

In  the  troposphere  clouds,  C02  and  H „0  control  the  neat  budget  of  the  atmos¬ 
phere  through  their  absorption  and  emission  of  infrared  radiation.  Other 
infrared  active  constituents  such  as  CH„,  N0o,  No0  and  CO  are  produced 
near  the  surface  of  the  earth  in  trace  amounts.  At  high  altitudes,  the  action 
of  solar  radiation  and  subsequent  chemical  reactions  produces  species  such 
as  Os,  NO,  NO+,  N20,  H202,  OH,  and  HOg. 

At  low  altitudes,  the  predominant  mechanism  leading  to  infrared  emission 
is  thermal  excitation  of  vibrational  and  rotational  energy  levels  of  the  ground 
electronic  state  of  these  minor  constituents.  At  high  altitudes,  however, 
other  activation  processes  must  be  considered. 

The  relative  importance  of  each  depends  on  the  molecule  involved  and  on  the 
altitude.  The  additional  mechanisms  known  to  be  of  importance  are  chemical 
formation  in  an  excited  state,  fluorescent  emission  due  to  excitation  by  visible 
or  ultraviolet  solar  irradiation,  resonant  excitation  by  solar  or  terrestrial 
infrared  radiation,  and  excitation  by  corpuscular  radiation  in  the  normal  day¬ 
time  ionosphere  or  during  aurorae.  Some  of  these  processes  continue  day 
and  night,  as  does  thermal  excitation,  while  others  are  restricted  to  daytime 
or  to  distrubed  atmospheric  conditions. 

In  addition  to  the  radiation  from  polar  molecules  in  excited  vibrational  levels 
of  the  ground  electronic  state,  transitions  between  e:-  cited  electronic  states 
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of  N2,  and  provide  a  sour  ce  of  infrared  emission.  Figure  23  shows 
the  position  of  infrared  vibration-rotation  bands  of  selected  molecular  species. 
The  long  lines  indicate  the  position  of  the  centers  of  strong  bands.  The  short 
lines  show  the  band  centers  of  weaker  vibration -rotation  bands  except  for  NO, 
where  they  indicate  the  position  of  transitions  between  highly  excited  electronic 
states.  Band  positions  were  obtained  from  Herzberg  (1945),  Wallace  (1962a, 
1962b)  and  Wray  (1969). 


B.  CHEMICAL  ABUNDANCES 

Very  little  experimental  data  is  available  on  the  distribution  with  altitude  of 
minor  constituent  molecular  abundances.  Only  ozone  and  nitric  oxide  con¬ 
centrations  have  been  determined  directly  at  altitudes  as  high  as  100  km. 
Hence  reliance  must  be  made  on  computational  estimates  which  are  only  as 
good  as  the  reaction  rates  used.  Several  important  rates  are  known  only  to 
within  an  order  of  magnitude,  while  Nicolet  (1965a)  points  out  that  really 
meaningful  computations  require  20  percent  accuracy. 

The  most  complete  computations  of  chemical  abundances  at  high  altitude  are 
those  of  Hunt  (1966)  and  Keneshea  (1967).  Accurate  computations  must  in¬ 
clude  the  effects  of  eddy  and  molecular  diffusion.  Unfortunately,  such  trans¬ 
port  mechanisms  have  as  yet  not  been  effectively  coupled  with  chemical 
computations  for  minor  atmospheric  constituents.  Computed  and  observed 
concentrations  may  differ  by  as  much  as  an  order  of  magnitude  as  is  the  case 
of  nitric  oxide  between  60  and  90  km. 


1.  Carbon  Dioxide 


The  most  abundant  minor  constituent  of  the  atmosphere  is  carbon  dioxide. 

Its  abundance  is  well  known  in  the  lower  atmosphere,  and  is  relatively  stable 
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Rotational  Spectrum 


at  about  320  parts  per  million  except  close  to  industrial  areas  where  the  burn¬ 
ing  of  fuel  increases  its  concentration  locally.  At  altitudes  above  about  70  km 

O 

it  is  dissociated  by  sunlight  in  the  spectral  range  1200  to  1750A.  Although 
Bates  and  Witherspoon  (1952)  gave  a  qualitative  discussion  of  the  chemistry 
of  CO 2  at  high  altitude  over  15  years  ago,  there  has  been  no  quantitative  study 
of  its  abundance  at,  high  altitudes.  Because  of  its  importance  as  a  source  of 
infrared  radiation,  knowledge  of  its  abundance  at  high  altitudes  is  required. 

As  part  of  the  current  effort,  CO2  photochemical  equilibrium  concentrations 
above  60  km  altitude  have  been  computed  using  solar  continuum  and  Lyman- 
Alpha  fluxes  from  Johnson  (1965).  The  most  important  reactions  considered 
are: 


CO2  +  hu  -* 

CO  +  O 

(5) 

CO  +  OH 

C02  +  H 

(6) 

CO  +  O  +  M 

-*  C02  +  M 

(7) 

The  reaction  rate  for  (6)  was  taken  from  Baulch  et  al.  (1968).  Reaction  (7) 
has  been  studied  by  many  experimental  groups,  with  widely  differing  results. 
The  rate  given  by  Barth  (1964)  was  considered  representative  and  was  used  in 
the  calculations. 

Because  of  the  shielding  of  CC^  from  the  continuum  flux  by  molecular  oxygen, 
Lyman -Alpha  radiation  at  1216A  accounts  for  most  of  the  photodissociation  of 
CC>2.  Even  so,  the  lifetime  of  C02  against  photodissociation  is  on  the  order 
of  a  year  and  atmospheric  diffusion  will  be  important.  On  this  basis  the 
photochemical  equilibrium  results  were  used,  neglecting  Lyman-Alpha 
radiation.  The  results  of  this  computation  are  given  in  Figure  24.  This 
abundance  distribution  should  be  considered  a  lower  limit  to  C02  abundance. 
Also  shown  as  a  dashed  line  is  the  CO 2  concentration  profile  given  by 
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Bortner  and  Kummler  (1968).  This  may  be  considered  an  upper  limit  to  the 
abundance  of  CQ9  with  altitude. 


2.  Water  Vapor 


Water  vapor  concentrations  vary  greatly  near  the  surface  of  the  earth. 

There  have  been  no  comprehensive  measurements  above  balloon  altitudes. 
Measurements  from  balloons  give  widely  varying  values,  and  it  is  probable 
that  water  vapor  carried  up  with  the  balloon  strongly  affects  measurements 
at  the  highest  accessible  altitudes.  Recent  stratospheric  observations  are 
those  reported  by  Mastenbrook  (1968).  In  a  series  of  51  flights,  observed 
volume  concentrations  in  the  lower  stratosphere  were  between  2  and  5  parts 
per  million  in  nearly  all  cases,  with  a  mean  of  about  3  parts  per  million. 
Extrapolation  of  such  measurements  to  higher  altitudes  is  uncertain,  but 
there  is  now  general  agreement  that  this  is  the  correct  order  of  magnitude  to 
use  in  chemical  computations  at  altitudes  up  to  100  km.  The  results  of 
Hunt  (1966)  are  adopted  here.  He  uses  a  similar  volume  mixing  ratio  of  5 
parts  per  million  as  the  initial  starting  point  of  his  chemical  computations. 
Hunt's  final  number  densities  of  HgO  are  shown  in  Figure  25. 

The  important  reactions  which  produce  HgO  at  high  altitudes  are: 


H  4  HOg  -  HgO  +  o 

(8) 

OH  +  HOg  HgO  +  Og 

(9) 

OH  L  OH  HgO  +  O 

(10) 

OH  +  HgOg  -  HgO  +  HOg 

(11) 

OH  +  Hg  HgO  +  H 

(12) 
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A  sixth  reaction  which  produces  water  vapor  is: 


O  +  CH4  -  H20  +  CH2  (13) 

Bates  and  Nicolet  (1965)  have  discussed  this  and  show  it  to  be  an  important 
source  in  the  stratosphere,  accounting  for  as  much  as  half  of  the  observed 
water  vapor  concentrations.  Hunt's  neglect  of  this  reaction  is  not  important 
if  the  FI20  produced  by  it  is  considered  as  part  of  the  initial  concentration. 


Those  reactions  which  destroy  H2Q  are: 


H20  +  hj  -  II  f  oil 

(14) 

H20  +  0(!D)  _  OH  +  OH 

(15) 

3.  Ozone 


1  he  other  molecules  treated  by  Hunt  which  are  considered  as  possible  sources 
of  infrared  emission  are  O^,  H202,  OH,  and  HOg.  Hunt's  abundances  for 
these  species  are  shown  in  Figures  26  through  28.  The  ozone  concentrations 
agree  wpII  with  measurements  below  80  km  and  will  not  be  appreciably  affected 
by  the  uncertainty  in  water  vapor  concentration.  The  concentrations  of  the 
other  species  must  be  considered  as  estimates  only,  although  the  OH  production 

rates  agree  well  with  the  observed  altitude  profiles  of  the  OH  infrared  band 
systems. 


Heestvedt  (1965)  has  done  a  similar  computation  for  several  latitudes  and  two 
values  of  the  water  vapor  mixing  ratio.  This  work  is  useful  as  a  guide  to 
possible  seasonal  and  latitudinal  variation  of  the  ozone  concentration,  but  is 
limited  to  giving  daytime  photochemical  equilibrium  values  of  concentrations. 
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The  abundance  of  H202  is  nearly  three  orders  of  magnitude  less  than  that  of 
Og •  The  vibration-rotation  bands  of  OH  are  at  wavelengths  shorter  than  5  /uni. 
Although  H02  is  present  at  80  km  in  about  the  same  abundance  as  O^,  there 
is  at  present  insufficient  data  available  to  make  meaningful  calculations. 

These  three  molecules  are  therefore  omitted  from  further  consideration  in 
the  current  study. 


4.  Nitric  Oxide 

The  daytime  concentration  of  nitric  oxide  is  shown  in  Figure  29.  This  plot 
is  the  result  of  combining  daytime  observations  by  Barth  (1966)  and  Pearce 
(1969)  with  computational  estimates  by  Nicolet  (1965b),  Ghosh  (1968)  and 
Saxena  (1968).  The  nighttime  concentration  should  not  be  appreciably  different 
at  altitudes  between  70  and  120  km  (Barth,  1965a  and  Saxena,  1968).  Below 
70  km,  nitric  oxide  is  destroyed  at  night  by  reaction  with  ozone.  Saxena 
(1969)  has  included  the  effects  of  diffusion  on  night  abundances  and  obtained 
a  lowering  of  concentration  above  120  km.  This  information  was  obtained 
too  late  to  include  in  the  computations. 

Before  Barth  (1964)  first  measured  nitric  oxide  concentrations  at  high  alti¬ 
tudes,  it  had  been  assumed  that  the  principal  reaction  forming  nitric  oxide  is 

N(4S)  +  02(V)  -  NO(2IIu)  +  0(3P)  (16) 

where  all  reactants  and  products  are  in  the  ground  state.  The  measured 
abundances  are  an  order  of  magnitude  greater  than  those  which  can  be  generated 
in  this  way  at  altitudes  below  90  km.  Subsequently,  Nicolet  (1965b)  proposed 
the  reaction 

°2  (2^g}  +  N2(lyg}  -  NO(2IIu)  +  NO+(1y+)  (17) 
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However,  the  measured  upper  limits  to  the  rate  make  this  reaction  unimportant 
below  100  km.  Hunten  and  McElroy  (1968)  have  proposed  the  reaction 

N(4S)  +  Q2  (!Ag)  -  NO  (2IIU>  4  0(3P)  (18) 

DelGreco  and  Kenneally  (1968)  have  reported  a  lower  limit  for  the  reaction 

- 14  3 

rate  to  be  4  x  10  cm  /sec.  Burt  and  Schiff  (1969),  however,  find  a  rate 

_  i  s  q 

constant  less  than  10  cm  /sec. 

The  lower  limit  to  the  rate  measured  by  DelGreco  and  Kenneally  (1968)  is 
about  25  times  too  small  to  produce  the  number  densities  found  by  Pearce 
(1969)  below  90  km.  Another  possible  reaction  is 

N(2D)  +  Q2  (3vg)  _  NO  (2IIu)  +  0(3P)  (19) 

Its  importance  depends  on  the  extent  to  which  dissociative  recombination  of 
NO  forms  atomic  nitrogen  in  the  D  state.  At  present  there  is  no  reaction 
known  with  certainty  to  be  capable  of  producing  observed  nitric  oxide  con¬ 
centrations  below  90  km.  However,  Geissler  and  Dickenson  (1968)  have 
examined  transport  phenomena  near  the  mesophere  and  conclude  that  nitric 
oxide  produced  above  85  or  90  km  can  be  carried  down  rapidly  enough  to 
explain  known  abundances  between  70  and  85  km.  Observed  abundances  are 
adopted  as  the  basis  for  the  computations. 

Nicolet  (1965a)  estimates  daytime  abundances  of  N02  to  fall  rapidly  with  in¬ 
creasing  altitude.  At  70  km  N09  concentrations  are  10  3  that  of  NO  and  at 
-4  ^ 

100  km  are  10  that  of  NO.  N02  concentrations  at  night  approach  the  NO 
daytime  concentration  below  75  km  due  to  the  reaction 

NO  +  03  -  N02  +  02  (  20) 

but  maintain  the  daytime  value  above  75  km. 
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5.  Nitrous  Oxide 


Nitrous  oxide  (N20)  concentrations  have  been  computed  by  Keneshea  (1967) 
and  are  shown -in  Figure  30.  For  noon  near  the  equator,  the  principal 
reaction  leading  to  formation  of  N20  is 

N2  +  O  +  M  _  N20  +  M  *21) 

The  principal  loss  reactions  are 

N20  +  hy  -»  N2  +  O  (22) 

N20  +  hy  -  N  +  NO  (23) 

N20  +  02(!A)  -  N2  +  °3  (24) 

There  is  no  significant  diurnal  variation. 


6.  Methane 


Methane  (CH^)  is  uniformly  mixed  in  the  troposphere  with  a  volume  mixing 
ratio  of  about  1.6  parts  per  million.  Figure  31  shows  concentrations  above 
50  km  on  the  assumption  that  the  mixing  ratio  remains  constant  in  the  upper 
atmosphere.  Cadle  and  Powers  (1966)  have  computed  the  rate  at  which  it  is 
destroyed  by  reaction  (19).  At  60  km  the  lifetime  may  be  as  short  as  a  day 
and  is  less  than  a  month.  Figure  31  is  thus  at  best  an  upper  limit  to  methane 
abundance. 
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ALTITU 


Figure  30.  Nitrous  Oxide  Concentration  Profile 
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7.  Other  Species 


Other  species  which  may  provide  a  source  of  infrared  radiation  are  oxides 

such  as  MgO,  CaO,  NaO  and  SiO  which  are  deposited  in  the  upper  atmosphere 

by  meteors.  Donohue  and  Meier  (1967)  have  discussed  the  chemistry  of  Na 

in  the  upper  atmosphere,  basing  their  arguments  on  the  observed  diurnal 

variation  of  the  sodium  dayglow.  They  conclude  that  influx  rates  are  such 

13  1 

that  sodium  in  all  forms  is  present  at  about  2  x  10  molecules /cm  above 

4  3 

50  km.  If  mixing  is  uniform,  the  number  density  would  be  5  x  10  /cm'  at 
90  km.  The  abundances  of  other  metallic  oxides  and  SiO  is  expected  to  be  of 
the  same  order  of  magnitude.  "Phis  is  three  orders  of  magnitude  below  the 
concentrations  of  NO  and  N90.  Therefore,  these  species  have  been  omitted 
from  further  consideration  in  the  current  study. 


C.  SUMMARY 

All  molecules  present  in  the  atmosphere  except  N^,  O 2  and  are  emitters 
of  infrared  radiation.  In  the  spectral  range  5  to  25  gm,  the  most  important 
at  altitudes  above  70  km  will  be  0O2,  U90,  Og,  NO,  NT20  and  possibly  CH4, 
by  virtue  of  their  abundance  relative  to  other  species. 

It  is  not  possible  at  present  to  give  exact  values  of  concentrations  as  a 
function  of  altitude.  Observational  data  exists  only  for  O^  and  NO.  For 
other  species  the  results  of  computations  must  be  used.  Such  computations 
are  only  as  good  as  the  input  data  and  physical  model  used.  In  the  case  of 
NO,  computations  and  observations  do  not  agree,  probably  because  transport 
processes  have  not  yet  been  reliably  included  with  the  chemical  computation 
scheme. 

A  rough  estimate  of  the  reliability  of  the  altitude-concentration  profiles  given 
in  this  section  is  as  follows.  C09  concentrations  are  bracketed  by  the  two 
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curves  shown  in  Figure  24.  Water  vapor  concentrations  may  be  a  factor  of 
two  lower  or  a  factor  of  five  higher  than  indicated  in  Figure  25.  Ozone  is 
strongly  controlled  by  the  photochemistry  and  equatorial  values  are  almost 
certainly  within  a  factor  of  two  of  the  values  of  Figure  26.  Nitric  oxide 
concentrations  below  120  km  will  be  within  a  factor  of  three  of  those  shown 
in  Figure  29  below  120  km,  with  expected  deviation  to  lower  and  not  higher 
values.  Above  120  km,  concentrations  will  lie  between  the  curves  for  day 
and  night.  In  view  of  the  known  discrepancy  between  computation  and  ob¬ 
servation  for  NO,  the  nitrous  oxide  curve  of  Figure  30  may  be  in  error  by  a 
factor  of  three.  No  estimate  is  possible  for  methane,  which  is  not  expected 
to  be  present  in  appreciable  quantities  above  60  km, 

In  view  of  the  uncertainties  in  atmospheric  abundances,  the  computational 
model  developed  in  this  study  and  described  in  the  next  section  is  expected 
to  be  more  accurate  than  present  knowledge  of  abundances  and  capable  of 
determining  abundances  to  greater  accuracy  when  used  to  analyze  experimental 
data  from  probes. 
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SECTION  V 

DEVELOPMENT  OF  COMPUTATIONAL  MODEL 


This  section  of  the  report  presents  the  mathematical  development  of  a 
realistic  physical  model  for  computing  infrared  radiances  from  gaseous 
molecules  in  the  upper  atmosphere.  The  abundances  of  important  molecules 
as  a  function  of  altitude  were  discussed  in  Section  IV,  and  these  data  will  be 
used  as  input  information  to  the  computational  model.  A  computer  program 
was  developed  during  the  study  to  perform  calculations  as  specified  by  the 
mathematics  of  the  model.  The  program  has  the  capability  of  performing 
exoatmospheric  calculations  of  limb  radiances  for  tangent  heights  from  60 
to  500  km  and  endoatmospheric  calculations  as  a  function  of  zenith  angle 
from  0°  (looking  up)  to  90°  (horizon  view)  and  for  any  reference  altitude 
within  the  atmosphere.  The  physical  approximations  assumed  in  the  model 
generally  hold  for  radiance  calculations  at  altitudes  above  70  km. 


A.  GENERAL  CONSIDERATIONS 

The  vibrational  population  of  molecules  determines  the  radiation  from  the 
molecules.  Several  factors  influence  the  vibrational  population;  collisional 
excitation,  radiative  excitation  and  fluorescent  excitation  processes,  all 
of  which  are  balanced  by  radiative  de -excitation  and  collisional  de -excitation. 

In  the  troposphere  and  stratosphere  collisional  excitation  of  vibrational 
levels  of  polar  molecules  is  sufficiently  rapid  to  compensate  for  radiative 
de -excitation  and  the  vibrational  population  levels  are  in  thermodynamic 
equilibrium  with  translation  and  rotation.  At  altitudes  above  70  or  80  km, 
however,  radiative  transitions  become  more  effective  than  collisions  in 
determining  vibrational  populations.  Absorption  of  radiation  becomes  in 
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many  cases  the  most  important  process  leading  to  infrared  radiation.  At 
wavelengths  longer  than  about  5  |im,  the  most  intense  source  of  radiation  is 
the  thermal  emission  from  either  carbon  dioxide,  water  vapor  and  ozone 
from  lower  levels  in  the  stratosphere,  troposphere  or  the  earth's  surface. 
During  the  day,  however,  solar  radiation  should  be  included  because 
fluorescent  processes  involving  absorption  of  visible  or  ultraviolet  radiation 
may  lead  to  appreciable  excitation  of  vibrational  levels.  For  molecules  with 
vibration-rotation  bands  at  wavelengths  shorter  than  6  or  7  pm,  solar  induced 
excitation  may  be  comparable  to  that  caused  by  absorption  of  surface  or 
atmospheric  radiation. 

Knowledge  of  vibrational  populations  of  molecular  species  in  the  upper 
atmosphere  is  essential  to  the  computation  of  band  and  spectral  radiances. 

A  major  effort  of  the  current  study  has  been  the  development  of  mathematical 
model  which  treats  in  combination  the  effects  of  collisional  excitation  and 
de -excitation  of  vibrational  levels  and  the  absorption  and  re -emission  of 
radiant  energy.  The  significant  new  result  of  this  study  is  a  simple  method 
of  treating  the  radiative  transfer  problem  for  a  single  molecular  vibration- 
rotation  band  with  Doppler  line  shape.  This  method  allows  the  treatment  of 
an  entire  band  as  if  it  were  a  single  line  and  greatly  simplifies  the  computa¬ 
tion  of  the  population  of  vibrational  states  of  a  molecvie  under  conditions 
where  there  is  a  lack  of  thermodynamic  equilibrium. 

The  molecular  species  considered  in  the  computational  model  include  CC>2, 
H20,  03,  NO,  N20  and  CH4.  Goody  (1964)  provides  a  good  introduction  to 
the  detailed  structure  and  infrared  radiative  properties  of  these  molecules. 

The  following  discussion  of  this  section  treats  the  mathematical  development 
of  the  computational  model,  including  the  physical  conditions  of  absorption 
and  re -emission  of  radiant  energy  when  absorption  can  be  neglected,  the 
treatment  of  radiative  transfer  in  the  presence  of  absorption,  collisional 
excitation  mechanisms,  and  the  details  of  determining  infrared  radiation 
from  the  molecules. 
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B.  BASIC  RADIATIVE  PROCESSES 


This  subsection  discusses  the  physics  of  the  interaction  of  light  with  an 
isolated  atom  or  molecule.  An  understanding  of  the  basic  processes  is  the 
first  step  in  any  computation  of  the  radiation  from  any  gas  or  mixture  of 
gases.  The  first  paragraph  introduces  the  problem  through  use  of  the 
Einstein  A  and  B  coefficients  and  relates  these  to  the  experimentally  deter¬ 
mined  oscillator  strength  and  band  strengths.  Succeeding  paragraphs  treat 
absorption  from  a  continuum  source,  a  spectrally  varying  source  and 
fluorescent  processes.  Later  subsections  will  discuss  the  modification 
necessary  when  a  large  number  of  molecules  are  present,  the  effects  of 
collisions  on  exciting  molecules  to  statesTrom  which  radiation  may  occur 
and  the  combination  of  radiative  and  collisional  processes  in  a  computational 
model  for  estimating  radiances  from  C09,  H20,  Og,  N20  and  CH4  at  high 
altitudes. 

When  treating  a  problem  in  which  the  number  of  photons  emitted  or  absorbed 
by  an  atom  or  molecule  must  be  considered  rather  than  total  energy  emitteo, 
it  is  natural  and  convenient  to  examine  the  physics  of  the  interaction  of 
radiation  with  matter  through  use  of  the  Einstein  A  and  B  coefficients.  In 
treating  a  radiative  problem  such  as  the  current  study  of  high  altitude 
infrared  radiances  which  involves  transfer  of  energy  between  molecules 
through  collisions,  this  approach  is  almost  a  necessity. 

The  Einstein  A^_^  coefficient  is  the  rate  at  which  an  atom  or  molecule  in  an 
excited  state  spontaneously  emits  radiation,  in  units  of  photons/ sec- 
molecule.  In  an  assemblange  of  molecules  all  in  a  state  which  can  radiate, 
the  average  rate  of  radiation  from  each  molecule  is  given  by  the  Einstein 
Au-£  Coefficient  for  the  Particular  radiative  transition  from  a  state  of  higher 
energy  (subscript  u)  to  one  of  lower  energy  (subscript  t). 
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In  addition  to  spontaneous  emission,  two  other  processes  occur.  These  are 

absorption,  characterized  by  the  Einstein  coefficient  for  absorption  B£->u, 

and  induced  emission,  characterized  by  the  Einstein  coefficient  for  induced 

emission  B  .  .  In  a  radiation  field  with  spectral  density  p,  the 
u-*i  vJt  u 

probability  that  an  atom  or  molecule  absorbs  a  quantum  of  radiation  in  unit 

time  is  B,  p  .  When  the  radiation  density  is  expressed  in  C.  G.S.  units 
*—  u  v. 

i  U  2 

ergs -sec /cm,  the  coefficient  B^u  has  units  cm  /ergs -sec.  Similarly,  an 

atom  or  molecule  in  a  radiation  field  with  spectral  density  Pviu  and  in  an 

excited  state  emits  radiation  at  the  rate  B  .  9  in  addition  to  that 

n  1  Viu 

characterized  by  the  coefficient  A  ^  The  Einstein  coefficients  are  proper¬ 
ties  of  an  atom  or  molecule.  If  one  coefficient  is  known,  the  others  may  be 
determined  by  the  relations 


=  8rrh  B  . 

u-i 


(25) 


and 

^i®i-u  ^u  ®u -*l 


(26) 


Thus,  it  is  possible  to  treat  the  problem  knowing  only  the  Einstein  Au_^ 
coefficient. 


Experimental  determinations  of  the  interaction  of  infrared  radiation  with 
gases  are  usually  determined  by  measurements  of  the  absorption  of  light. 
Experimental  results  may  be  expressed  as  integrated  absorption 


8tt 


fu 

Si 


(27) 
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or  band  strengths 


S .  =  S'  /pc 

l u  £u  ^ 


(28) 


In  the  visible  and  ultraviolet  regions  of  the  spectrum,  it  is  customary  to 
use  oscillator  strengths  to  characterize  the  interaction  of  radiation  and 
matter.  In  this  case,  the  defining  equation  may  be  taken  to  be 


Tr  e 
me 


N.f„ 
l  £-u 


(29) 


The  results  of  theoretical  computations  of  the  interaction  of  matter  with 
radiation  are  frequently  expressed  as  dipole  moments.  These  are  related 
to  the  Einstein  A  coefficient  by 


A 


u-  l 


64tt4 

3h\3 


(30) 


The  rate  at  which  a  single  atom  or  molecule  absorbs  and  reradiates  light 
from  a  parallel  beam  from  a  continuum  source  at  a  wavelength  Xq  corres¬ 
ponding  to  a  resonance  transition  is  (Mitchell  and  Zemansky,  1961): 

2 

G  =  10  4  — I  f  X  2  (photons /sec) 

2  o  o  o 
me 

(31) 

=  8.852  x  10  *  I  f  X  3  (photons /sec) 

o  o  o 

2 

If  the  incident  flux  is  expressed  in  photons /cm  -sec -pm  when  the  incident 

2 

light  intensity  is  expressed  in  watts /cm  -  pm, 


-24  2 

G  =  3.  72  x  10  *  I  SX  (watts) 


(32) 


Both  of  these  equations  are  necessary  because,  in  the  details  of  the  model, 
photons  are  counted,  while  final  results  are  presented  as  spectral  power. 
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Wnen  the  source  is  an  extended  black  body  surface  with  a  spectral  radiance 
N.  ,  the  G  factor  is  obtained  by  integration  over  the  solid  angle  subtended  by 
the  source.  For  radiation  from  a  plane  parallel  surface  of  infinite  extent, 

-9^  2 

G  =  1.  169  x  10  NXS  X  (33) 

One  other  purely  radiative  mechanism  is  of  importance  in  upper  atmosphere 
radiance  computations.  This  is  fluorescence,  and  must  be  considered  in 
treating  radiation  from  CO,,,  NO,  NgO  and  CH^.  A  striking  example  occurs 
for  the  weak  10.  4  pm  COQ  band.  Absorption  of  solar  radiation  by  the  C09 
4.  3  pm  band  maintains  the  vibration  mode  at  a  vibrational  temperature  of 
about  290°K.  During  the  day,  this  leads  to  an  enhancement  of  emission  in  the 
10.  42  pm  band  by  as  much  as  two  orders  of  magnitude  when  viewed  along  the 
earth's  limb  from  outside  the  atmosphere. 


E<luorescence  occurs  when  the  upper  state  of  the  resonance  transition  can 
reradiate  into  states  other  than  the  initial  one.  The  relative  probability  of 
a  radiative  transition  to  a  lower  state  n  from  the  upper  state  u  when  there 
are  k  possible  lower  states  is: 


P 


u-n 


u-*n 


k 


J  =  1 


(34) 


where  the  A  are  the  Einstein  A  coefficients  for  the  spontaneous  emission 
u-J 

probabilities  for  a  transition  from  the  upper  state  to  the  jth  lower  state. 

In  the  case  of  molecular  resonance -fluorescent  scattering  at  low  temperatures, 
only  the  lowest  vibrational  level  (v"  =  0)  in  the  ground  electronic  state  is  sig¬ 
nificantly  populated.  Absorption  of  a  photon  excites  the  molecule  to  a  vibrational 
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level  v'  in  the  same  or  a  higher  electronic  state.  Radiation  is  possible  from 
this  state  to  the  initial  level  (v"  =  0)  or  to  other  vibrational  levels  (v"  /  0)  in 
the  ground  electronic  state.  The  rate  of  resonance  fluorescent  scattering  in  a 
given  electronic  vibrational  transition  may  be  obtained  by  combining 
Equations  (31)  and  (34)  and  introducing  appropriate  subscripts: 


G  ,  „  =  8.  853  x  10"17  I  .  f  ,  \  , 
v'v  v'o  v'o  v'o 


A 


v’v" 


E 


A  |  M 
v'v 


(3  5) 


Here  the  subscript  v'v"  refers  to  a  transition  from  an  upper  state  with  vibra¬ 
tional  level  v'  to  a  lower  state  with  vibrational  level  v".  The  subscript  v'o 
refers  to  absorption  of  a  photon  by  a  molecule  in  the  lowest  vibrational  level 
of  the  ground  state,  exciting  the  molecule  to  a  vibrational  level  v'  in  the 
same  or  a  higher  electronic  state. 


C.  RADIATIVE  TRANSFER  FOR  LARGE  OPTICAL  PATHS 

/ 

In  developing  a  realistic  model  of  the  interaction  of  radiation  with  matter, 
it  is  necessary  to  take  into  account  the  details  of  absorption.  Line  shape 
must  be  considered  because  it  determines  the  amount  of  radiation  escaping 
from  an  optically  thick  layer  of  gas.  Thus,  for  large  optical  thickness 
(greater  than  a  magnitude  of  10),  the  total  emission  from  an  isolated  spectral 
line  with  Lorentz  shape  is  proportional  to  the  square  root  of  the  logarithm  of 
the  optical  thickness.  In  addition,  absorption  of  radiation  emitted  at  one 
place  in  the  gas  by  a  molecule  at  another  place  cannot  be  computed  without 
knowledge  of  the  line  shape. 
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-At  high  altitudes  the  profile  of  the  spectral  absorption  curve  of  an  individual 
line  in  a  vibration -rotation  band  is  predominantly  that  due  to  velocity  or 
Doppler  broadening.  Kuhn  and  London  (1969)  have  investigated  the  heat 
budget  of  the  atmosphere  between  30  and  110  km.  They  report  that  radiative 
transfer  calculations  require  the  use  of  the  Voigt  profV  (including  both 
collisional  and  Doppler  broadening)  at  altitudes  be. .wet  ;  20  and  70  km.  but 
that  the  Doppler  profile  is  adequate  at  higher  alti-LuGos  J  e  Doppler  shape 
being  the  largest  contributor  to  line  broadening  move  ,  m  for  CO.,,  H,(J 
and  03<  The  computational  model  therefore  assumes  i  Doppler  profile  for 
simplicity. 

In  addition  to  the  use  of  Doppler  shape  for  a  single  line,  several  other  as¬ 
sumptions  are  made  about  the  radiative  transfer  within  lines  of  a  band. 

These  are  the  following: 

(1)  Rotational  levels  within  a  vibrational  level  remain  in  equilibrium 
with  the  translational  temperature.  This  is  valid  up  to  at  least 
150  km  according  to  Goody  (1964). 

(2)  There  is  complete  redistribution  of  frequency  within  a  single 
line  following  absorption  of  radiation. 

(3)  There  is  complete  redistribution  of  absorbed  energy  among  all 
rotational  levels. 

(4)  The  variation  of  line  strength  with  rotational  level  is  taken  to  be 
that  given  by  the  Hb'nl-London  formulae,  neglecting  wavelength 
variation  across  a  band. 

(5)  There  is  no  overlap  of  adjacent  lines. 

(6)  Temperature  variations  at  different  levels  do  not  affect  the 
radiative  transfer.  This  last  assumption  is  probably  the 
largest  source  of  error  in  computing  transmission. 


73 


The  chief  result  is  to  make  computed  band  radiances  too  small  when  deter¬ 
mining  escape  of  radiation  from  a  warm  layer  overlaid  by  a  cooler  layer. 

It  will  be  shown  in  Section  VI  A(H90  limb  radiances)  that  this  approximation 

(Lt 

yields  an  insignificant  error  in  computed  radiances  for  tangent  heights 
greater  than  75  km. 


The  above  assumptions  lead  to  the  follDwing  treatment,  which  begins  with 
consideration  of  a  single  line  and  from  this  develops  the  method  used  for 
treating  radiative  transfer  in  an  entire  band.  The  spectral  absorption 
coefficient  of  a  Doppler  broadened  line  is  {Mitchell  and  Zemansky,  1961): 


k(v)  = 


rre 

me 


N/o 


(36) 


wnere  the  Doppler  half  width  at  half  maximum  is 


/  2  RT  Ini 
*  Me2 


v  (sec  *) 
o 


The  amount  of  light  absorbed  by  the  line  per  unit  frequency  interval  is 


I(v)  A(v)  dv  =  I(v)  (1  -  T(v)  dv 

=  I(v)  (1  -  e"k(v))  dv 


(37) 


and  the  total  light  absorbed  is 


00 

J  I(v)  A(v)  dv 
o 


J  I(v>  (1  -  T(v)  )  dv 
o 


=  J  I(v)  (1  -  exp 
o 


(38) 
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where  the  integration  is  over  all  frequencies  and 


K  =  —  A/ 
0  AvdV 


In  2  n  e 


tt  me  ^i'o 


The  line  will  absorb  only  near  the  center  frequency  v  and  I(v)  may  be 
replaced  by  I(vq).  Making  the  substitution 


v  -  v 
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The  corresponding  expression  on  a  wavelength  basis  is 


(39) 
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(40) 
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T1  e  quantity  in  parentheses  in  Equation  (40)  is  simply  a  factor  by  which  the 
product  of  G  and  must  be  multiplied  to  account  for  the  effects  of  total 
optical  thickness  in  computing  the  total  absorption  of  radiation  by  a 
column  containing  molecules /cm^  in  the  lower  state.  This  quantity  is 
called  S  by  Mitchell  and  Zemansky  (1961)  who  present  tabulations  of  its 
value.  Penner  (1959)  gives  an  asymtotic  expansion  for  large  values  of  N  . 
The  astrophysical  literature  frequently  uses  the  notation 


2  ot 


L(Ko)  =  Ko  S(«Q)  =—  J 

TT  O 


(41) 


The  factor  by  which  G  must  be  multiplied  to  obtain  the  true  rate  of  absorption 
at  an  optical  depth  may  be  found  by  differentiating  Equation  (41)  with 
respect  to  (Ivanov  and  Shcherbakov,  1965): 


M,  (  K  ) 
l  o 


d  L(KQ) 
d  K 


2  ao  2 

p  -w 

—  J  e 

TT  O 


K  p“w 


dw 


(42) 


The  factor  by  which  G  must  be  multiplied  to  obtain  the  actual  absorption  from 
a  black  body  source  is  (Ivanov  and  Shcherbakov,  1965) 

2tt  rr/2  K 

IW  =  J  dcp  J  cos  9  sin  9  — — —  d0  (43) 


o 


o 


cos  0 


When  the  source  of  radiation  is  a  Doppler  profile  exactly  that  of  the  absorbing 
line,  the  differential  absorption  at  optical  thickness  is  given  by 
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e-w 


dw 


(44) 


M2(«0  ) 


rr  o 


2 


and  G  must  be  multiplied  by  this  factor.  Also,  if  the  source  is  an  infinite 
plane  layer  (not  a  surface)  radiating  with  a  Doppler  profile,  G  must  be 
multiplied  by 

tt/2  /  K  \ 

N2(}V  =  J  dcP  J  M2{— —  )  sin  q  dp  (4b) 

o  o  \  cos  0/ 

The  functions  of  Equations  (41)  through  (45)  were  evaluated  by  numerical 
integration  and  polynomial  approximations  obtained  for  the  entire  range  of 
variation  of  K0  from  zero  to  infinity,  using  the  known  asymtotic  forms  as  a 
guide  at  large  values. 

If  the  line  strengths  for  each  line  of  a  vibration-rotation  band  of  a  linear 
molecule  such  as  CC>2,  NO,  and  NgO  are  determined,  and  the  optical  thickness 
of  the  strongest  line  in  the  band  is  taken  as  a  parameter,  functions  analagous 
to  those  in  Equations  (41)  to  (45)  may  be  computed  numerically.  It  is  noted 
that  temperature  affects  the  Doppler  width  in  the  calculations.  However, 
results  indicate  that  a  temperature  variation  of  up  to  a  factor  of  50  results 
in  less  than  a  10  percent  variation  in  the  values  of  these  functions.  Therefore, 
only  two  sets  of  computations  are  required  for  linear  molecules,  one  for 
perpendicular  bands  and  one  for  parallel  bands.  These  two  sets  are  adequate 
even  in  cases  where  nuclear  spin  effects  cause  alternation  of  intensities  of 
rotational  lines. 

The  method  of  extending  these  functions  to  an  entire  band  is  illustrated  for 
the  function  (KQ)  of  Equation  (42)  for  a  parallel  -  ^Y  transition. 
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The  rotational  energy  levels  of  a  linear  molecule  in  a  ^  state  are  given  to  a 
first  approximation  by 


E(J)  =  J(J  +  1)  By  he  (46) 

The  relative  populatic  of  a  single  rotational  level  is  (2J  +  1) 
exp(-E(J)/kt)/Q^  where  is  the  rotational  partition  function 


Q  =  T  (2 J  +  1)  exp(-E(J)/kt)  (47) 

R  J  =0 

The  band  for  a  * 7  transition  contains  two  sets  of  lines  or  branches,  one 

corresponding  to  a  change  in  J  between  the  upper  and  lower  state  by  +1  (R 
branch),  the  other  to  a  change  in  J  of  -1  (P  branch).  The  relative  strength 
of  each  line  is 

R 

S  =  (J"  +  1)  (48) 

J 


Sp  =  J"  (49) 

J 

At  a  given  temperature  the  (fictitious)  J  value  of  the  strongest  line  is  computed 
by  setting  the  derivative  with  respect  to  J  of  (J"  +  1)  exp  (-E(J)/kt)  equal  to 
zero 


J"max  V~  -1  (5f» 

■tJ  v  f  1 

v  he 

An  optical  thickness  Kq  is  assigned  to  the  strongest  line  and  the  optical 

thicknesses  for  all  lines  in  the  band  are  computed.  The  product  K  M1  (K) 

is  computed  for  each  line  and  summed  for  all  lines  for  which  the  product  is 
-14 

more  than  10  that  of  the  strongest  line.  This  sum  is  divided  by  the  sum  of 
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r 


optical  thicknesses  to  obtain  the  band  function 


M 


1 


B 


v  v 

J  I  =  P,  R 


T  V 

J  I  =  P,  R 


(51) 


Similar  band  functions  corresponding  to  Equations  (41),  (43)  to  (45)  are 
indicated  by  using  a  superscript  B  with  the  symbol  for  the  basic  function. 
Polynomial  approximations  have  been  obtained  for  all  the  functions  and  are 
used  in  the  numerical  calculations. 


Time  has  not  permitted  a  similar  treatment  of  the  more  complex  bands  of 
^2^'  ^3  C H ^ .  A  proper  treatment  will  require  the  use  of  a  complete 
set  of  line  strengths  such  as  that  being  developed  by  Dr.  McClatchey  at 
AFCRL.  In  treating  these  non-linear  molecules  for  the  current  investigation, 
a  fictitious  value  of  Bv  has  been  estimated  which  leads  to  the  effective  use  of 
about  500  lines  per  band  instead  of  the  thousands  actually  present. 


D.  INCLUSION  OF  COLLISIONAL  EXCITATION  IN  THE  MODEL 

In  computing  radiances  at  high  altitudes,  it  must  be  recognized  that 
vibrational  levels  are  not  populated  at  thermal  equilibrium  values. 

Radiation  depopulates  vibrational  states  rapidly  and  collisional  excitation 
is  not  efficient  enough  to  maintain  a  Boltzmann  distribution  of  vibrational 
states.  Collisional  excitation  and  de- excitation  rates  are  required  to 
compute  the  balance  of  opposing  mechanisms  which  determine  actual 
vibrational  populations,,  This  subsection  discusses  the  addition  of  collisional 
mechanisms  to  the  radiative  processes  described  previously  and  illustrates 
the  use  of  the  model  for  water  vapor  radiance  computations. 

/' 

i 


79 


Three  collisional  processes  are  considered  which  affect  the  vibrational 
state  of  a  molecules.  These  are:  (1)  translational- vibrational  (T-V) 
interactions  in  which  translational  energy  is  transformed  into  vibrational 
energy  or  vice-versa  during  a  collision,  (2)  intermolecular  vibrational 
exchange  (V-V)  in  which  a  molecule  of  one  species  exchanges  a  quantum  of 
vibrational  energy  with  a  molecule  of  another  species  and  (3)  intramolecular 
vibrational  exchange  in  polyatomic  molecules  in  which  a  collision  transfers 
the  energy  in  one  mode  of  vibration  into  another  mode  within  the  same 
molecule. 


When  a  molecule  in  the  ground  state  collides  with  another  molecule  there  is 
a  finite  probability  Py^  that  a  part  of  the  relative  translational  energy  will 
be  transformed  to  vibrational  energy  if  the  translational  energy  is  equal  to 
or  exceeds  the  vibrational  energy.  Similarly,  if  the  molecule  is  in  an 
excited  vibrational  state,  there  is  a  finite  probability  P^y  that  the 
vibrational  energy  will  be  converted  to  relative  translational  energy  during 
the  collision.  For  a  Boltzmann  distribution  of  translational  energies  at  a 
given  temperature. 


P 


TV 


^u 

-  PyT -  exp 


(52) 


Similar  expressions  hold  for  intermolecular  and  intramolecular  V-V  trans¬ 
fer. 


The  rate  at.  lich  a  vibrational  level  is  excited  (krpy)  or  de-excited  (kyy) 
is  given  by  the  product  of  the  collision  frequency  Z  and  the  transition 
probability.  Experimental  determinations  of  collisional  excitation  rates  may 
be  reported  as  probabilities,  excitation  rates  or  relaxation  times.  The 
latter  are  related  to  excitation  and  de-excitation  rates  by  (Herzfeld  and 
Litovitz,  1959) 
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1 


( 53 ) 


T 


k  4-  k 

VT  TV 


Excitation  rates  depend  on  temperature  and  number  density  through  the  strong 
temperature  variation  of  the  transition  probability  and  the  square  root 
temperature  and  linear  density  variation  of  the  collision  frequency. 

All  parts  of  the  model  are  now  ready  for  use.  The  computational  method 
may  be  illustrated  by  application  to  a  system  in  which  only  one  excited 
vibrational  level  is  considered,  say  water  vapor.  This  level  is  assumed  to 
be  excited  and  de-excited  by  both  radiative  and  collisional  processes. 

O 

At  a  given  altitude,  there  will  be  N  molecules /cm  in  the  ground  state  and 
Nu  molecules/cm  of  nitrogen  and  oxygen  which  are  the  dominant  collision 
partners  in  collisional  excitation  and  de-excitation  processes. 

In  unit  volume,  molecules  in  the  ground  state  are  excited  to  the  upper  state 

by 

(1) 

(2) 

(3) 

(4) 

(5) 


Collisional  T-V  excitation  at  a  rate  N  k 

Xj  1  V 

Collisional  V-V  excitation  at  a  rate  N.k,.,, 

I  M  V 


Absorption  of  solar  flux  at  a  rate  N.I  G'  M  B  (X  ) 

I  o  1  '  s 

Absorption  of  lower  atmosphere  radiance  at  a  rate 
Ni\,TG'NlB  «r> 

Absorption  of  radiance  emitted  at  other  levels  of  the 
atmosphere  at  a  rate  N  f  N,  „  G  '  NB  (K  )  dZ 

£  K,  U  Z  U 
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Molecules  in  the  upper  state  are  de-excited  by 


(1)  Collisional  V-T  de-excitation  at  a  rate  N  k.rrp 

u  VT 

(2)  Collisional  V-V  de-excitation  at  a  rate  N  k.r,, 

U  V  IVl 

(3)  Spontaneous  radiative  de-excitation  at  a  rate 
For  each  level,  there  is  a  differential  equation 


dN  R  R 

— —  =  N£  (kTV  +  kMV  +  IQ  C’  Mx  (Ks)  +  N  G'Ni  (kt) 
dt 

+  r  \  D  G'  N2B  <V  dZ)  ‘  N£  (kVT  +  kVM  +  AuM  '  0  <54) 


The  equations  are  solved  iteratively,  with  starting  conditions  being  those 
appropriate  to  thermodynamic  equilibrium*  For  molecules  in  which  more 
than  one  vibrational  level  is  considered,  an  equation  must  be  added  at  each 
altitude  level. 

After  the  vibrational  populations  have  been  determined  at  each  level,  both 
exospheric  and  endospheric  radiances  are  computed,  taking  into  account 
the  necessary  geometry. 

A  brief  comment  on  treatment  of  nitrogen  vibrational  temperatures  is  in 
order  at  this  point.  At  altitudes  above  100  km,  it  is  possible  for  the  nitrogen 
vibrational  levels  to  be  populated  above  the  thermal  equilibrium  condition, 
either  by  electron  impact  or  as  a  result  of  collisional  quenching  of  atomic 
oxygen  in  the  excited  *D  electronic  state.  Walker  (1968)  has  discussed  the 
excitation  of  N2  resulting  from  the  quenching  of  0(*D).  His  analysis  of 
Langmuir  probe  measurements  of  E-region  electron  temperatures  concludes 
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that  at  110  km  and  above,  the  vibrational  temperature  of  N2  is  of  the  order 
of  3  000°Kor  more  during  the  daytime.  When  the  program  was  set  up  to 
include  such  high  vibrational  temperatures,  it  was  found  that  the  nitrogen 
transferred  more  energy  to  the  lower  density  C02  concentration  model 
near  110  km  than  could  be  produced  by  quenching  of  O^D).  This  suggests 
that  C02  may  play  an  important  role  in  relaxing  the  N9  vibrational 

Lt 

temperature  at  altitudes  up  to  120  km. 

The  program  contains  the  option  of  using  either  ambient  translational 
temperature  or  an  arbitrary  temperature  profile  for  the  N9  vibrational 

Li 

temperature. 

The  following  subsection  describes  the  molecular  data  used  for  each  of 
the  species  CC>2,  H2Q,  C>3,  MO,  N20  and  CH4. 
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L.  APPLICATION  OF  MODEL  TO  INDIVIDUAL  MOLECULAR  SPECIES 

I  his  subsection  discusses  briefly  the  data  which  were  used  in  the  band 
radiance  computations  reported  in  Section  \  I.  Table  1  summarizes  values 
of  band  strengths  and  Table  2  summarizes  collisional  excitation  rates.  It 
must  be  recognized  that  the  latter  are  not  as  well  known  as  the  former.  In 
some  cases  the  exact  path  through  which  vibrational  energy  is  transferred 
is  not  known.  The  vibrational  energy  level  scheme  used  for  each  molecule 
is  treated  below. 


1 .  Carbon  Dioxide 

Figure  32  shows  the  vibrational  energy  le\  els  that  have  been  included  in  the 
computations,  together  with  the  first  excited  vibrational  energy  level  of 
molecular  nitrogen.  Seven  levels  in  addition  to  the  ground  state  have  been 
considered,  and  nine  radiative  transitions,  eight  of  which  lie  between  10  and 
17  g m.  Hie  scheme  is  a  good  first  approximation  to  the  real  physical  proc¬ 
esses.  Levels  formed  by  combinations  of  vibrational  modes  have  been  ne¬ 
glected  but  are  considered  of  less  importance.  Dashed  lines  indicate  radia¬ 
tive  transitions  and  solid  lines  indicate  paths  by  which  collisions  transfer 
vibrational  energy  or  transform  translational  to  vibrational  energy.  The  high 
band  strength  for  C02  at  4.  3  microns  shown  in  Table  1  is  expected  because 
the  frequency  is  high  and  the  vibration  is  a  fundamental. 


Collisional  excitation  rates  are  taken  from  Bitter  man  (1969).  The  Y-T  rates 
for  the  transition  between  the  ground  state  and  level  2  were  applied  to  the 
other  transitions  in  the  y2  vibrational  series,  making  the  necessary  modifi¬ 
cations  due  to  the  differing  statistical  weights  of  the  levels.  Band  strengths 
for  the  seven  bands  near  15.0  gm  were  taken  from  Drayson  and  Young  (1966). 
Those  for  the  10.4  gm  and  4.3  gm  bands  were  taken  from  Burch  et  al.  (1962). 
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Table  1.  Band  Strengths  of  Infrared  Vibration -Rotation  Bands 


Molecule 

Band  Tenter 
wavelength,  |im 

— 

Band  Strength 
(ATM-cnri  )  -1 

Temperature 

References 

C°2 

4.3 

2500 

;  299 

Burch  et  al.  (1962) 

10.  4 

0.023 

299 

Burch  et  al.  (1962) 

13.9 

5.  0 

300 

Drayson  and  Young 
(1  967) 

15.  0 

194 

300 

Drayson  and  Young 
(1967) 

15.  0 

15.  0 

300 

Drayson  and  5  oung 
(1967) 

15.  0 

0.  85 

300 

Drayson  and  Young 
(1967) 

15.4 

1.0 

300 

Drayson  and  Young 
(1967) 

16.  2 

4.27 

300 

Drayson  and  Y'oung 
(1967) 

16.  7 

0.  14 

300 

Drayson  and  Young 
(1967) 

h2o 

6.  3 

300 

273 

Ludwig  et  al.  (196  5) 

°3 

9.  1 

10.  6 

273 

Clough  and  Kneizys 
(1967) 

9.  6 

455 

273 

NO 

5.3 

122 

273 

n2o 

4.  5 

1850 

303 

Burch  and  Williams 
(1962) 

7.  8 

12 

303 

Burch  and  Williams 
(1962) 

17.  0 

33 

303 

Burch  and  Williams 
(1962) 

ch4 

3.3 

320 

2  73 

Burch  and  Williams 
(1962) 

5.  8 

0.  01 

300 

Estimate 

6.  6 

2.  5 

2  73 

Burch  and  Williams 
(1962) 

7.  7 

185  | 

273 

Burch  and  Williams 
(1  962) 

NOTE:  Nitric  acid  has  absorption  bands  at  7.  5  and  11.  2  microns;  however, 
band  strengths  were  not  available  at  the  writing  of  this  report. 
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Table  2.  Collisional  Excitation  Rates 


A.  Vibrational- Translational  Energy  Exchange  Probabilities 


Molecule 

p 

VT 

Temperature 

Reference 

co2 

3. 6  x  10‘ 6 

300 

Taylor  and  Bitterman 
(1969) 

h2o 

1  x  10"  ° 

300 

Taylor  and  Bitterman 
(1969) 

°3 

1  x  10'4 

300 

Estimate 

NO 

4  x  10" 8 

292 

Estimate 

n2o 

7  x  10"  6 

297 

Yardley  (1  968) 

ch4 

9  x  10"4 

296 

Yardley  and  Moore 
(1968) 

B.  Vibrational- Vibrational  Energy  Exchange  Probabilities 


Molecule 

p 

mv 

Temperature 

Reference 

co9 

2. 3  x  10~3 

300 

Taylor  and  Bitterman 

(1969) 

H„0 

1  X  10‘5 

300 

Taylor  and  Bitterman 

L. 

(1969) 

°3 

-4 

5x10 

300 

Estimate 

NO 

4x10"' 

292 

Basco  et  al.  (1961) 

n2° 

8. 7  x  10"4 

297 

Yardley  (1  968) 

ch4 

2.  3  x  10"3 

296 

Yardley  and  Moore 

(1968) 
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2.  Water  Vapor 


Because  the  water  vapor  molecule  is  extremely  efficient  in  relaxing  its 
vibrational  excitation,  vibrational  states  other  than  that  which  gives  rise  to  the 
6.3  jum  band  have  been  neglected.  As  for  C02,  vibrational  relaxation  and 
exchange  data  were  taken  from  the  curves  of  Taylor  and  Bitterman  (1969). 

I'he  band  strength  used  was  taken  from  Ludwig  et  al.  (1965). 


3 .  Ozone 

No  experimental  data  on  vibrational  energy  transfer  in  ozone  were  found. 
Values  for  vibrational -translational  and  vibrational-vibrational  transfer 
intermediate  between  those  for  methane  and  water  vapor  have  been  used. 
The  two  bands  at  9.  1  pin  and  9.6  pirn  were  treated  as  a  single  band,  using 
the  total  strengths  given  by  Clough  ana  Kneizys  (1967). 


4.  Nitric  Oxide 

The  vibrational-translational  energy  exchange  rate  for  nitric  oxide  in  nitro¬ 
gen  has  not  been  determined  accurately  because  the  rate  for  vibrational- 
vibrational  exchange  is  much  faster.  The  experimental  determination  of  the 
latter  by  Basco  et  al.  (1961,  1962)  has  been  used  and  the  former  taken  to  be 
one -tenth  that  value.  The  band  strength  for  the  fundamental  was  taken  from 
Abels  and  Shaw  (1966)  and  that  for  the  first  overtone  from  Schruin  and  Ellis 
(1966). 

Eluorescent  excitation  following  absorption  of  ultraviolet  solar  radiation  by 
the  gamma  band  system  at  wavelengths  below  2270A  has  been  included  as 
follows.  Neglecting  variation  of  solar  flux  and  transition  probability  with 
wavelength,  the  Einstein  A  coefficients  are  proportional  to  the  Eranck-Condon 
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factors  qv,vM.  1  he  relative  population  of  each  vibrational  level  v'  of  the 
upper  electronic  state  is  given  by  the  Franck-Condon  factor  q  ,  (assuming 
absorption  only  from  the  lowest  vibrational  level  v"  =  0  of  the  ground 
electronic  state).  For  each  value  of  v',  the  sum  v  t  v"  q  ,,  gives  the  mean 
vibrational  excitation  of  the  molecule  after  emission  of  visible  or  ultraviolet 
radiation.  the  double  sum  ^  qv,o  v"  qv, v»»)  gives  the  mean  vibrational 
level  of  the  molecule  following  all  stiort  wavelength  excitations.  Using 
Franck-Condon  factors  from  Nicholls  (1964),  the  mean  vibrational  level  of 
NO  after  excitation  is  3.  85. 

khe  gamma  band  absorption  was  approximated  by  a  single  gand  with  an 
oscillator  strength  the  average  of  that  of  the  v1  =  0  to  v1  =  4  bands. 

Chemiluminescence  from  the  reaction 

N  f  02  NO  +  O 


with  both  reactants  in  the  ground  state  was  included  using  the  production  rate 
curve  of  Ghosh  (1968).  The  average  vibrational  level  during  formation  is 
assumed  to  be  v"  =  4,  as  suggested  by  Dalgarno  (1963). 


5.  Nitrous  Oxide 

Figure  33  shows  the  energy  level  scheme  adopted  for  N20.  As  for  C02, 
vibrational  excitation  processes  are  shown  by  solid  lines  and  radiative  proc-' 
esses  by  dotted  lines.  The  upper  multiplets  in  the  u2  were  treated  as  single 
levels.  Einstein  A  coefficients  were  assigned  according  to  the  statistical 
weights  of  the  upper  levels  (Equation  34),  and  all  radiations  assumed  to  be  at 
the  fundamental  wavelength.  Collisional  excitation  rates  were  obtained  from 
Yardley  (1968),  who  states  that  the  experimental  data  do  not  permit  a  good 
determination  of  the  rate  of  population  of  the  ^  vibrational  level.  Level  6 
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Figure  33. 


n2o 


Nitrous  Oxide  Vibrational  Levels  and  Collisional 
and  Radiative  Transitions  Included  in  Infrared 
Radiance  Model 
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has  been  coupled  to  level  3  by  an  interaction  with  the  same  strength  as  that 
for  exchange  of  vibrational  energy  between  nitrogen  and  levels  5  and  7. 


6.  Methane 

Figure  34  illustrates  the  energy  level  scheme  considered  for  methane.  Again, 
experiment  does  not  give  a  complete  picture  of  vibrational  exchange  rates. 
Interactions  with  the  infrared  inactive  Uj  fundamental  have  been  neglected. 

The  collisional  excitation  data  used  were  obtained  from  Yardley  and  Moore 
(1968)  and  band  strengths  from  Burch  and  Williams  (1962). 


F.  SUMMARY 

This  section  of  the  report  has  described  the  development  of  a  mathe¬ 
matical  model  which  can  be  used  to  calculate  infrared  band  radiances  in  the 
upper  atmosphere  for  the  molecules  C02,  H20,  O^,  NO,  N20  and  CH^. 

The  physical  assumptions  and  the  expected  accuracy  of  the  computational 
model  are  commensurate  with  uncertainties  in  the  knowledge  of  physical 
properties  of  the  upper  atmosphere. 

The  model  incorporates  both  collisional  and  radiative  excitation  mechanisms 
and  collisional  and  radiative  de -excitation  mechanisms  when  computing  the 
populations  of  vibrational  states.  Collisional  mechanisms  are  considered  in 
as  great  a  detail  as  available  experimental  results  permit.  In  addition, 
fluoi  °scent,  and  chemiluminescent  mechanisms  are  included  as  processes 
which  populate  the  vibrational  states.  For  those  altitude  regions  of  the 
atmosphere  where  radiati\e  transfer  can  be  important  (about  70  to  100  km), 
the  model  utilizes  a  Doppler  line  shape  when  calculating  the  transmission 
properties  of  the  atmosphere.  The  assumption  that  rotational  lines  in  a 
band  do  not  overlap  makes  it  possible  to  treat  the  radiative  transfer  problem 
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6.5yM 


v  2 


Figure  34. 


Methane  Vibrational  Levels  and  Collisional 
and  Radiative  Transitions  Included  in  Model 
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for  an  entire  band  of  a  linear  molecule  as  simply  as  if  it  were  a  single  line, 
t  inally,  the  detailed  computational  scheme  for  each  molecule  is  outlined  for 
future  reference. 
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SECTION  VI 

HI  SI  LTS  I  "HO  VI  COMPUTATIONAL  MODEL 


This  section  of  the  report  presents  the  results  of  computations  of 
high  altitude  infrared  radiances  using  the  model  described  in  the  previous 
section.  The  calculations  are  based  on  the  six  species  COg,  H2^  °3J 
N  O  and  CH.  which  are  expected  to  contribute  most  significantly  to  high- 
altitude  radiation  in  the  5  to  25  jum  spectral  region.  The  best  estimate 
abundance  profiles  of  Section  IV  have  been  used  with  temperature  and 
pressure  profiles  from  the  U.S.  Standard  Atmosphere,  1962.  Upwelling 
terrestrial  and  lower  atmosphere  intensities  at  the  60  km  level  were 
obtained  from  the  low  altitude  radiance  program  described  in  Section  II. 

Solar  flux  values  were  obtained  from  Johnson  (1964). 

Band  radiance  values  for1  noon  and  night  are  presented  for  all  six  species  for 
exoatmospheric  viewing,  and  endoatmospheric  calculations  are  illustrated  by 
band  radiances  for  the  nitric  oxide  5.  3  gm  band.  Spectral  radiances  irom 
the  water  vapor  rotational  band  are  given  for  exoatmospheric  viewing.  It  is 
emphasized  that  these  results  are  estimates  based  on  uncertain  chemical 
abundances.  Certainly  further  investigations,  including  radiance  measure¬ 
ments  by  probes  and/or  satellites  will  be  required  to  determine  temporal 
and  spatial  variations  of  radiances  and  information  concerning  chemical 
abundances  of  atmospheric  species. 


A.  LIMB  RADIANCES 

1 .  Carbon  Dioxide 

At  present  the  computer  program  calculates  the  vibrational  populations  of  the 
seven  vibrational  levels  of  C09  shown  in  Figure  3  2  and  calculates  band  radiances 
for  nine  bands,  seven  of  which  lie  near  15  jum. 
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It  is  of  interest  to  examine  the  effects  of  each  of  the  exciting  mechanisms 
Vt hich  c  ontribute  to  maintaining  the  populations  of  the  various  levels  and 
hence,  the  resulting  radiation  field.  The  results  in  Figure  35  illustrate 
these  effects  and  show  the  sum  of  the  limb  radiances  of  the  seven  bands  near 
15  fim  for  the  highest  expected  C()9  concentration.  Curves  are  drawn  for 
ioui  conditions;  local  thermodynamic  equilibrium  (LTL),  collisional  excitation 
only  (but  including  radiative  transport  of  radiation  arising  from  collisional 
excitation),  nighttime  illumination  from  the  lower  atmosphere,  and  equatorial 
noon  illumination.  The  CO^  v9  fundamental  vibration  mode  which  radiates  at 
15.0  /im  is  populated  efficiently  by  collisions,  and  the  first  excited  vibrational 
level  remains  populated  within  a  factor  of  two  ol  the  thermal  equilibrium  value 
up  to  about  70  km.  1  his  is  not  apparent  from  the  curves  in  1*  igure  35  when 
comparing  LTE  and  the  collisional  excitation  curve  because  the  optical  depth 
is  large  along  a  limb  view  path  and  much  of  the  radiation  arises  from  altitudes 
above  the  tangent  height.  A  large  fraction  of  the  radiation  from  lower* 
altitudes  is  absorbed  by  C09  at  high  altitudes  before  escaping  to  space.  When 
the  radiation  from  the  earth  and  lower  atmosphere  is  added,  there  is  only  a 
small  increase  in  band  radiance,  noticeable  on  the  scale  of  the  drawing  only 
at  tangent  heights  between  85  and  105  km.  The  noon  curve  shows  a  larger 
increase  because  some  of  the  solar  energy  absorbed  by  the  strong  Vn  funda¬ 
mental  mode  at  4.  3  g m  is  converted  to  excitation  of  vibrational  levels  which 
radiate  in  the  15  gm  spectral  region. 


Figure  36  shows  the  same  data  for  the  band  at  10.4  gm.  The  large  increase 
in  radiance  (one  to  three  orders  of  magnitude)  above  70  km  at  noon  is  due  to 
fluorescent  emission  following  absorption  of  solar  energy  at  4.3  gm. 

Figures  37  and  38  give  LTE,  noon  and  night  radiances  for  the  lower  limit  to 
C02  abundance  (see  figure  24).  For  the  15  gm  band  complex  and  the  10.42 
band,  respectively,  there  is  little  difference  in  radiances  for  tangent  heights 
below  90  km.  The  radiance  magnitudes  associated  with  the  low  abundances 
are  one  to  two  orders  of  magnitude  less  than  those  shown  for  the  high  abundances 
profiles  (Figures  35  and  36;  for  tangent  heights  greater  than  100  km. 


BAND  RADIANCE  NAX  (WATTS/CM  STER) 


60  70  80  90  100  110  120 

TANGENT  HEIGHT  h  (KM) 


Figure  3  5.  Limb  Radiance  of  Carbon  Dioxide  15  urn.  Bands 
Upper  Limit  Concentration  Profile 
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Figure  38.  Limb  Radiance  of  Carbon  Dioxide  10.  4  Mm  Band, 
Lower  Limit  Concentration  Profiles 
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2,  Water  Vapor 


Figure  39  shows  a  comparison  of  the  6.3^  H20  limb  view  radiances  computed 
for  LTE,  noon  and  night  conditions.  Below  7  5  km,  collisional  processes 
maintain  the  vibrational  population  of  the  v2  niode  first  vibrational  level  near 
equilibrium.  Between  75  and  95  km,  the  absorption  of  radiation  from  lower 
levels  of  the  atmosphere  maintains  the  vibrational  population  above  that  which 
holds  at  LTE,  even  at  night.  At  noon,  absorption  of  solar  radiation  increases 
the  radiance  by  about  a  factor  of  two  over  that  of  night  at  all  altitudes  above 
75  km. 

It  is  noted  that  the  computational  model  used  to  calculate  band  radiances  from 
non-linear  molecules  represents  a  first  approximation  to  the  absorption 
problem  because  the  number  of  lines  is  underestimated  (see  Section  V).  The 
lower  dashed  line  at  the  left  of  the  illustration  is  for  LTE  conditions  with  an 
extremely  large  over -estimation  of  the  effects  due  to  absorption.  As  can  be 
seen,  the  radiance  magnitudes  are  depressed  by  an  order  of  magnitude  for 
tangent  heights  of  60  to  65  km.  However,  there  is  little  difference  between 
the  two  LTE  curves  for  tangent  heights  greater  than  75  km.  it  is  concluded 
from  this  comparison  that  the  treatment  of  absorption  by  non-linear  molecules 
in  the  computational  model  will  result  in  significant  errors  only  at  the  lower 
tangent  heights  where  the  physical  approach  is  marginal  anyway.  At  tangent 
heights  greater  than  70  km,  the  treatment  of  absorption  is  reasonably  accurate 
and  the  errors  introduced  into  the  calculations  are  of  little  significance  to  the 
overall  magnitudes. 


3 .  Ozone 

Figure  40  presents  band  radiances  for  the  9.6  jum  band  for  noon  and  pre¬ 
dawn  conditions.  The  results  indicate  that  there  is  less  than  an  order  of 
magnitude  difference  in  radiance  magnitudes  at  all  heights  due  to  diurnal 
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Figure  40.  Limb  Radiance  of  Ozone  9.6f(mBand 


variations.  Ozone  reaches  its  maximum  concentration  just  before  dawn  and 
it  might  be  expected  that  solar  illumination  would  suddenly  enhance  the 
radiance  field.  However,  the  nighttime  ozone  distribution  has  appreciable 
optical  thickness  and  radiative  transfer  from  the  lower  atmosphere  and 
within  the  ozone  distribution  are  important.  Therefore,  solar  illumination 
of  the  pre-dawn  ozone  maxima  with  noon  flux  values  increases  limb  radiances 
by  no  more  than  about  50  percent  at  any  one  altitude.  It  can  be  predicted  on 
this  basis  that  the  effect  of  photodissociation  as  the  sun  rises  will  be  such 
that  limb  radiance  values  will  relax  toward  the  noon  values  with  a  time  scale 
on  the  order  of  minutes  during  sunrise. 


4.  Nitric  Oxide 

During  the  day,  and  at  altitudes  above  90  km,  nitric  oxide  is  expected  to 
radiate  more  strongly  in  the  5  to  25  g m  region  than  any  of  the  other  molecules 
studied.  The  chief  source  of  this  radiant  energy  is  chemiluminescence.  The 
results  of  limb  radiance  computations  for  the  NO  5.3  gm  fundamental  band 
are  shown  in  Figure  41.  The  contribution  from  each  excitation  mechanism 
is  shown  as  was  done  in  the  case  for'  CO^  (see  Figure  35).  The  solid  curves 
give  computed  radiances  based  on  the  assumption  that  the  nitrogen  vibrational 
temperature  is  equal  to  the  translational  temperature.  Asa  comparison,  the 
dashed  curves  show  radiances  computed  using  the  translational  temperatures 
up  to  100  km,  a  linear  temperature  variation  between  100  and  1  10  km  and  a 
3000°K  temperature  above  110  km.  The  effects  due  to  vibrationally  excited 
nitrogen  are  apparent  between  80  and  120  km.  The  results  show  that  the 
radiation  field  is  the  strongest  mechanism  populating  the  vibrational  levels 
during  the  night  and  partial  noon  at  other  heights.  However,  during  full  day, 
chemiluminescence  is  the  major  source  of  radiation,  and  the  uncertainty  in 
the  estimate  of  the  temperature  of  vibrationally  excited  nitrogen  becomes 
irrelevant. 
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Figure  41.  Limb  Radiance  of  Nitric  Oxide  5.3  Mm  Band 
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It  lb  emphasized  that  nitric  oxide  radiances  vary  more  between  day  and  night 
than  do  those  of  any  other  species  considered  in  this  study. 


5.  Nitrous  Oxide 

C  omputed  limb  radiance  proiiles  for  the  N00  7.8  and  17.0  /tm  bands  are 
shown  in  figures  42  and  43.  Radiation  from  the  lower  atmosphere  is  the 
dominant  source  of  excitation  for  the  7.  8  jum  band  at  all  altitudes  and  is  the 
most  important  exciting  mechanism  for  the  17.0  n m  band  above  80  km.  There 
is  practically  no  day-night  variation  in  radiances. 

6.  Methane 

Although  methane  is  not  expected  to  be  present  in  appreciable  amounts  cbove 
the  stratosphere  because  of  the  rapid  r  eaction  with  atomic  oxygen,  limb 
radiance  computations  have  been  performed  using  the  tropospheric  mixing 
ratio  as  an  upper  limit  estimate.  Figures  44,  15  and  46  present  the  results 
of  limb  radiance  calculations  for  the  CH4  5.  8  nm,  6.5  fira  and  7.  7  gm  bands. 
Absorption  of  radiation  is  the  chief  excitation  mechanism  for  all  of  these 
bands,  and  collisional  excitation  plays  a  negligible  role. 


B.  ENDOATMOSPHERIC  RADIANCES 

The  computational  model  developed  in  this  study  has  the  capability  of  per¬ 
forming  radiance  calculations  as  a  function  of  zenith  angle  from  an  assumed 
reference  platform  within  the  atmosphere.  Endoatmospheric  calculations 
are  useful  in  simulating  radiances  that  would  be  measured  by  rocket  borne 
radiometers  during  ascent  and  descent  portions  of  flight  programs.  Endo¬ 
atmospheric  radiances  were  calculated  for  all  the  bands  discussed  in 
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Section  A  above.  Only  one  example  is  included  in  the  current  report  since 
this  study  has  emphasized  limb  radiances  that  would  be  measured  from  exo- 
atmospherie  platforms. 

Figure  47  illustrates  the  variation  in  NO  band  radiances  as  a  function  of  zenith 
angle  and  reference  altitudes  of  100  km,  200  km  and  300  km.  The  effect  of 
increased  optical  paths  with  increasing  zenith  angle  is  apparent  in  the  in¬ 
creased  radiance  magnitudes  at  all  altitudes.  The  curvature  of  the  radiance 
variations  at  100  km  is  less  than  those  indicated  at  200  and  300  km.  This  is 
due  to  the  fact  that  the  transmission  properties  for  the  horizon  view  at  100  km 
are  becoming  optically  thick,  and  absorption  is  significant  and  radiative 
transfer  is  becoming  important.  The  calculations  for  the  higher  altitudes 
are  optically  thin.  Thus,  the  radiance  variations  with  zenith  angle  are 
similar  at  200  and  300  km,  a  function  principally  of  the  geometry  of  view. 


C.  SPECTRAL  RADIANCES 

The  final  goal  of  this  study  of  upper  atmospheric  infrared  radiance  is  the 
prediction  of  spectral  radiance  values  at  wavelengths  between  5  and  25  /.on 
for  all  species  and  for  both  exoatmospheric  and  endoatmospheric  viewing. 

As  an  example,  Figure  48  shows  spectral  limb  radiances  for  the  water  vapor 
rotational  band  at  wavelengths  between  10  and  25  grn  for  tangent,  heights  of 
60  km,  80  km  and  100  km.  These  rad.ances  were  computed  by  a  program 
which  is  included  as  a  subroutine  in  the  final  radiance  model  computer  pro¬ 
gram.  The  basic  assumption  used  in  the  computation  is  that  rotational  levels 
are  in  thermodynamic  equilibrium  and  that  emission  and  absorption  may  be 
computed  using  Equations  42  and  44.  The  necessary  dipole  moments  were 
obtained  from  the  dipole  matrix  element  values  tabulated  by  Gates  et  al. 
(1964)  for  wavelengths  longer  than  12  gm  and  from  values  of  absorption  co¬ 
efficients  supplied  by  Dr.  R.  McClatchey  of  AFCRL  for  shorter  wavelengths. 
About  250  of  the  strongest  line^  m  this  spectral  region  are  included  in  the 
computational  results.  The  spectral  radiance  for  all  species  is  given  in 
Appendix  F. 
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Figure  47.  Endoatmospheric  Radiances  of  Nitric  Oxide  5.3  Mm  Band 
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D.  SUMMARY  OF  RESULTS  AND  COMPARISON  WITH  OBSERVATIONS 


A  summary  comparison  of  limb  radiance  profiles,  presented  in  Section  A, 
is  shown  in  Figure  49  for  equatorial  noon  conditions  and  for  important 
molecular  bands  considered  in  this  study.  The  results  presented  here  allow 
a  comparison  of  the  relative  magnitudes  of  different  bands  for  tangent  heights 
from  60  to  about  120  km.  It  is  noted  again  that  the  physical  assumptions 
made  in  the  computational  model  are  valid  above  the  70  km  tangent  height. 

By  far  the  dominant  source  of  radiance  at  all  altitudes  is  the  5.  3  gm  NO 
band.  The  15  g m  C02  complex,  representing  the  total  contribution  of  seven 
bands  between  13.9  and  16.7  gm  appears  to  be  the  second  largest  radiance 
source  at  heights  between  70  and  7  5  km,  and  at  heights  above  90  km  up  to 
where  the  molecule  is  dissociated.  CH^,  and  N90  have  about  the  same 
order  of  magnitude  in  the  70  to  100  km  altitude  range.  Above  100  km,  NgO 
is  about  an  order  of  magnitude  larger  than  the  radiances  from  Cll^  and  Og, 
as  well  as  that  from  the  10.4  gm  CC>2  band.  The  results  indicate  that 
radiances  from  the  6.  3  gm  C00  are  the  smallest  of  all  the  molecules  con- 
sidered  for  heights  above  about  75  km.  Furthermore,  the  radiance  magnitudes 
from  the  H20  rotational  spectra  in  Figure  48  indicate  the  same  order  of 
magnitude  as  those  for  the  6.3  gm  band.  Therefore,  1190  radiances  for  the 
abundances  assumed  in  the  upper  atmosphere  are  least  important  of  all 
species  considered. 

The  most  extensive  experimental  observations  which  correspond  to  the  exo- 
atmospheric  limb  radiances  presented  above  are  those  of  Markov  (1969). 

This  paper,  which  summarizes  several  articles  which  have  appeared  in  the 
Russian  literature,  reports  radiances  in  the  3  to  8  gm  spectral  region  which 
far  exceed  any  estimates  from  the  computational  model  of  the  current  study. 
Even  at  LTE  (Figure  41),  nitric  oxide  cannot  explain  Markov's  radiances, 
which  vary  from  10  to  7  x  10  watts /cm  (about  10  to  4  x  10  watts /cm  - 
ster)  in  layered  structure  at  altitudes  between  about  150  and  500  km. 
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Figure  49.  Summary  Comparison  of  Limb  Radiance  Profiles  for  Noon 
Conditions  and  Important  Molecular  Bands  Considered  in 
This  Study 
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The  most  plausible  partial  explanation  of  Markov's  results  invokes  collisional 
transfer  of  nitrogen  vibrational  excitation  to  carbon  monoxide  formed  by  photo¬ 
dissociation  of  CO 2  at  nigh  altitudes.  Excitation  of  CO  in  the  4.  7  g m  band  is 
about  200  times  as  efficient  as  excitation  of  NO  by  collision  (Walker,  196  8; 
Vasco  et.  al.,  1961).  Further,  there  may  be  up  to  10  times  as  much  CO  as  NO 

at  120  km.  This  could  result  in  limb  band  radiances  from  CO  on  the  order 
-5  2 

of  10  watts /cm  -ster  at  altitudes  near  120  km  which  is  about  equal  to  the 

lower  values  of  Markov's  observations.  CO  band  radiances,  however,  cannot 
explain  layered  structure  at  altitudes  above  120  km.  Balloon  observations  of 
spectral  radiances  performed  by  Hampson  and  discussed  by  Houghton  (1964) 
further  substantiate  the  hypothesis  that  CO  is  an  important  radiance  source 
in  the  upper  atmosphere. 

It  is  certainly  worthwhile  in  the  future  to  examine  CO  in  greater  detail  and 
to  estimate  its  contribution  to  the  infrared  radiant  background  of  the  earth. 
Other  molecules  such  as  N^O  and  OH  may  also  be  significant  radiance 
sources  at  wavelengths  below  5  /u m. 
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SECTION  VII 

INFRARED  RADIATION  FROM  PARTICULATE  MATERIAL 


The  results  presented  in  Section  \  I  indicate  the  order  of  magnitude  and 
variability  of  limb  radiances  from  gaseous  sources  in  the  upper  atmosphere. 
Limb  radiances  can  vary  by  five  orders  of  magnitude  for  different  wave¬ 
lengths  in  the  5  to  25  pm  spectral  region.  It  is  important  in  the  overall 
problem  of  determining  the  natural  infrared  radiant  background  of  the 
earth  to  estimate  the  spectral  radiances  arising  from  non -gaseous  sources 
such  as  particulate  material.  These  may  include  naturally  observed  phenomena 
such  as  clouds  and  dust  layers  as  well  as  infrared  radiation  from  dust  in  inter¬ 
planetary  space.  A  comparison  of  radiances  from  gaseous  sources  with  those 
from  particulate  sources  will  give  answers  to  questions  such  as  (1)  which 
source  is  dominant?,  (2)  over  what  regions  of  the  spectrum  is  each  source 
important?,  and  (3)  at  what  altitudes  is  a  given  source  important?  Answers 
to  these  questions  will  help  to  complete  the  overall  picture  of  the  earth's 
infrared  radiant  background. 

It  is  important  in  the  following  discussion  to  keep  in  mind  the  general  order 

of  magnitude  of  limb  radiances  from  gaseous  sources  in  the  80  km  to  100  km 

tangent  height  range.  Summarizing  the  magnitudes  presented  in  Section  VI, 

-7  2 

the  5.3  p m  NO  band  has  a  magnitude  of  about  10  watts/cm  -ster.  The 

2 

orders  of  magnitude  from  other  bands  in  units  of  watts/ cm w -ster  and  in- 

7 

creasing  wavelengths  are:  6.  3  pm  H^O  band  —  10  ;  6.  55  pm  CIl^  band  — 

6  x  10“2 * * * 6 7  to  5  x  10“12;  7.  66  pm  CH4  band  -  10_8  to  10~10 *;  7.  8  pm  N£)  band  - 

3x10  8  to  3  x  10  9;  9.  6  pm  Oq  band  —  2  x  10  8  to  10  9;  10. 4  pm  C09  band  — 

-0  -10  J  -13  Z 

2  x  10  1  to  2  x  10  ;  10-12  /  m  H90  rotational  band  -  10  ;  14-16  pm  CO^ 

band  -  3  x  10_B  to  5  x  10~9;  17  pm  N£0  band  -  3  x  10_f  to  2  x  10“J;  and  the 

20-25  pm  II^O  rotational  band  —  5  x  10  ^  to  5  x  10  It  is  noted  that 

regions  of  the  spectrum  having  low  spectral  radiances  are  those  in  the 
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- 1 3  2 

atmospheric  window  at  1 1  g m  (~  10  watts /cm“-ster-gm)  and  the  rotational 

- 1 0  2 

water  vapor  band  from  about  18  g m  to  25  g m  (~  10  watts /cm  -ster-gm). 
These  spectral  regions  should  be  kept  in  mind  when  considering  the  radiances 
from  particulate  material  sources. 


A.  GENERAL  CONSIDERATIONS 

The  infrared  background  radiation  between  5  and  25  g m,  as  observed  at  or  near 
the  earth's  limb  from  particulate  material,  can  be  produced  by  three 
mechanisms: 

(a)  Rayleigh  scattering  of  infrared  radiation  (particles  <  1  gm) 

(b)  Mie  scattering  of  infrared  radiation  (particles  >  1  gm) 

(c)  Thermal  emissions  from  particles. 

Infrared  radiation  from  interplanetary  space  is  considered  in  the  same  con¬ 
text  as  radiation  from  clouds  and  dust  layers  since  its  origin  is  from  par¬ 
ticulate  material  and  may  be  of  significance  in  the  absence  of  sources 
imbedded  in  the  earth's  atmosphere. 

The  two  contributions  due  to  scattering  above  require  an  outside  source  of 
energy.  The  important  sources  of  irradiation  are  the  earth-atmosphere 
system  and  the  sun.  The  magnitude  of  irradiance  from  the  earth  is  dominant 
in  the  upper  atmosphere  for  wavelengths  greater  than  about  6  gm.  On  the 
other  hand,  solar  irradiance  is  dominant  in  interplanetary  space  at  distances 
more  than  about  200,  000  km  from  the  earth. 

The  magnitude  of  infrared  radiation  from  particulate  material  depends  on  the 
wavelength  of  light,  the  number  of  particules,  the  size  distribution  of 
particles,  their  physical  characteristics  and  composition,  and,  in  the  case 
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of  thermal  emission,  the  temperature  of  the  particles.  I  he  first  step  in 
nr  holing  the  radiation  from  particulate  material  in  the  upper  atmosphere  and 
interplanetary  space  has  been  to  review  the  present  state  of  knowledge  of 
particle  number  concentration,  size  distribution  and  composition  (if  available), 
as  functions  of  altitude  above  the  surface  of  the  earth  from  about  20  km  on 
upwards.  Sixteen  references  were  reviewed  which  include:  Blandari  et  al. 
(1968);  Hodge  and  Wright  (1968);  Far  low  (1968);  Kossler  and  dePary  (1966); 
Kiaser  (1968);  f  riend  (1966);  Yewkirk  and  Eddy  (1964);  Heard  (1959); 

Divar  (1959);  Cadle  (1966);  .lunge  (1961);  Rosen  (1964);  Wright  and  Hodge 
(1968);  and  COSPAR  (1968).  This  literature  search  is  condensed  in  graphical 
j.orm  of  Figure  50  (as  far  as  particle  concentration  is  concerned)  which 
lepicts  the  particle  number  concentration  per  end3  as  a  function  of  altitude. 
Particle  sizes  covered  in  this  survey  range  from  about  0.  1  jum  to  about  10  pm 
in  diameter,  with  a  predominance  of  particles  of  the  order  of  1  pm  in  diameter. 
The  horizontal  lines  represent  concentration  ranges  either  measured  or 
estimated  by  particular  authors;  vertical  lines  represent  such  determinations 
over  a  given  range  of  altitudes;  and  the  data  points  represent  specific  measure¬ 
ments  or  estimates  at  a  given  altitude.  Arrows,  indicate  a  possible  continuation 
of  the  particular  line. 

figure  50  illustrates  the  enormous  range  of  uncertainties  and  discrepancies 
of  the  available  information  in  this  field,  even  at  altitudes  accessible  by  re¬ 
search  balloons.  If  any  significant  conclusion  can  be  reached  at  this  juncture, 
it  is  that  above  about  50  km,  or  even  lower,  the  dust  particle  concentration  falls 

Q 

off  very  rapidly  as  a  function  of  altitude  from  typical  levels  of  0.  1  to  1  cm 

at  30  km  to  less  than  10  cm  °,  at  about  100  km.  Beyond  this  altitude  the 

concentration  continues  to  decrease  to  a  level  about  which  there  is  at  best 

-  8  -3 

only  a  fair  agreement,  ranging  from  a  maximum  of  10  cm  to  a  minimum 

“.  1 1  —  3  16  3 

of  10  cm  (one  estimate  indicating  a  concentration  of  10  D  cm  ’’  appears 

_  C|  _3 

unreasonably  low)  with  a  most  probable  value  around  10  1  cm  '  .  If  the 
presently  preferred  model  of  a  heliocentric  zodiacal  dust  cloud  is  correct, 
its  contribution  to  the  infrared  background  radiation  can  be  significant.  Even 
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though  the  above-mentioned  number  concentration  of  particles  is  ver\  low  in 
interplanetary  space  at  one  astronomical  unit  from  the  sun,  the  path  lengths 
involved,  especially  along  the  ecliptic  plane,  are  considerable;  and  thus,  the 
integrated  contribution  of  radiation  can  be  significant. 

In  addition  to  the  relatively  invariant  or  stable  dust  concentrations  near  the 
earth  and  in  interplanetary  space,  there  exist  transient  and  localized  par¬ 
ticulate  concentrations  such  as  nacreous  clouds,  noctilucont  clouds, 
leuchtstreifen'  and  artificial  injections  mainly  due  to  the  combustion  oi 
solid  rocket  fuel  propellants. 

Nacreous  clouds  typically  appear  at  altitudes  below  30  km  and  will  not  be 
considered  for  the  present  study  since  their  efie<  t  is  overshadowed  by  gaseous 
emission  of  CO.-,,  and  Il9()  in  the  lower  atmosphere.  Insufficient  quantitative 
information  precludes  the  inclusion  of  either  the  leuchtstreifen  (luminous 
bands  at  about  120  to  130  km)  or  of  artifically  produced  noctilucent  clouds. 
Naturally  occurring  noctilucent  clouds  will  be  treated  in  detail  with  respect 
to  their  contribution  to  the  natural  infrared  radiant  background  of  the  earth. 


B.  RADIATION  FROM  THE  STABLE  DUS  1  CON  1  ENT  IN  THE  EARTH'S 
ATMOSPHERE 


As  pointed  out  above  (see  Figure  50)  particle  number  concentration  appears 
to  drop  off  very  rapidly  above  30  or  40  km  such  that  no  significant  con¬ 
tribution  of  infrared  radiation  can  be  expected  from  the  stable  dust  content 
above  those  altitudes.  Since  only  a  few  reliable  measurements  have  been 
performed  in  the  altitude  region  between  the  balloon  ceiling  and  low  orbit 
levels,  it  is  fair  to  assume  that  stable  particle  concentrations  in  that  region 
should  be  of  the  order  of  10  J  to  10"L  cm"1  ,  several  orders  of  magnitude 
below  the  stratospheric  values.  Preliminary  calculations  using  the  above 
concentrations  and  typical  associated  path  lengths  show  levels  of  infr;  red 
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radiance  well  below  those  obtained  from  other  radiating  mechanisms  and 
sources  even  in  10.3  gm  "window"  discussed  in  Section  III  of  this  report 
(see  subsequent  section  on  noctilucent  clouds). 

C.  RADIATION  FROM  INTERPLANETARY  DUST 

The  presence  of  dust  particles  outside  the  earth's  atmosphere  has  long  been 
surmised  on  the  basis  of  the  observation  of  Zodiacal  light.  The  characteristics 
of  this  phenomenon  has  led  investigators  to  attribute  it  to  the  presence  of  the 
so-called  Solar  E  corona  or  dust  corona.  'he  particles  associated  with  the 
dust  corona  appear  to  form  a  very  large  lenticular ly  shaped  cloud  around  the 
plane  of  the  ecliptic.  It  true  extent  is  unknown  but  its  outer  limits  may  extend 
as  far  as  the  orbit  of  Jupiter.  The  main  area  of  uncertainty  is  associated 
with  the  actual  number  density  of  the  particles  within  this  distributed  cloud. 
Estimates  have  ranged  from  10  ^  cm  ^  to  10  ^  cm  L  at  1  AU  from  the  sun, 
depending  on  whether  a  geocentric  or  a  heliocentric  zodiacal  dust  model  was 
preferred.  The  prevailing  trend  is  to  accept  the  latter  of  these  two  models 
with  typical  particle  number  concentrations  of  10  to  10  ^  cm  at  1  Al  . 
Measurements  of  infrared  radiation  from  zodiacal  dust  have  been  restricted 
to  solar  elongations  (angle  of  observation  with  respect  to  the  direction  of  the 
sun)  of  a  few  degrees,  and  at  wavelengths  shorter  than  the  band  of  interest 
of  5  to  25  /im,  e.g.,  2.  2  gm  (Peterson,  1967;  MacQueen,  1968).  Several 
authors  have  developed  models  of  infrared  emission  from  interplanetary  dust 
clouds,  among  them  are  Peterson  (1963)  and  Kiaser  (1968).  Kaiser's  model 
incot  porates  the  latest  experimental  observations  at  2.  2  gm  in  his  emission 
model  using  a  more  rigorous  treatment  than  Peterson.  Kaiser  utilizes  a 
composite  "dirty  quartz"  particle  model  and  extends  liis  calculations  to  13 
gm.  Again,  only  small  solar  elongation  angles  were  considered  (up  to  2.  5  ). 

Figure  51  shows  the  calculated  spectral  infrared  radiance  from  the  inter¬ 
planetary  dust  replotted  from  Kaiser's  work  (solid  lines).  These  curves 
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1  10  100 


ANGULAR  ELONGATION  FROM  SUN-DEGREES 


Figure  51.  Spectral  Radiance  from  Interplanetary  Dust 
for  Different  Elongation  Angles  from  the  Sun 
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have  been  extrapolated  (dashed  lines)  on  the  basis  of  the  general  trends  of 
Kaiser's  calculations.  These  extrapolated  curves  most  probably  furnish  an 
upper  limit  of  the  expected  infrared  radiance  from  interplanetary  dust  if  the 
assumptions  of  particle  number  concentration  and  size  distribution,  as  well 
as  composition  are  all  valid.  Figure  51  also  depicts  the  results  of  prior 
calculations  performed  by  Peterson  (1963),  for  5  gm  emission  and  scattering. 
The  scattering  contribution  falls  off  vary  rapidly  for  wavelengths  longer  than 
5  gm  and  can  thus  be  neglected  as  compared  with  the  expected  emission  levels. 

A  more  rigorous  treatment  of  the  emission  from  interplanetary  dust  in  the 
infrared  for  elongations  beyond  2.  5  was  not  attempted  at  this  point  in  view 
of  practical  time  constraints  and  based  on  the  confidence  which  the  overall 
trends  of  the  emission  characteristic  provide  to  a  simple  extrapolation  as 
shown  in  Figure  51.  The  validity  of  the  particular  model  of  the  spectral  charac¬ 
teristic  of  this  emission  is  predicated  on  a  possible  experimental  confirmation 
of  this  band  structure  associated  with  the  optical  characteristics  of  the 
material  assumed  to  predominate  in  the  form  of  interplanetary  dust  particles 
(Kaiser,  1968). 

The  results  in  Figure  51  indicate  that  the  magnitude  of  radiance  from  inter  plane 
tary  dust  decreases  with  increasing  angle  to  the  sun  and  depends  on  the  com¬ 
position  and  reststrahlen  bands  of  the  dust  material.  Spectral  radiance 
magnitudes  on  the  order  of  10  1  to  10  ^  watts/cm~ -ster-gm  appear  to  be 
possible  for  viewing  angles  of  about  50  degress  to  the  sun.  A  probable 
magnitude  for  angles  greater  than  90  degrees  is  between  10  and  10 

9 

watts /cnW -ster-gm.  These  results  are  in  general  agreement  with  recent 
calculations  made  by  Walker  (1969);  however,  they  do  not  agree  with  cal¬ 
culations  reported  by  Mikirob  (1966). 

It  is  noted  that  these  magnitudes  are  greater  than  limb  radiance  magnitudes 
at  11  gm  mentioned  earlier,  and  are  of  the  same  order  or  larger  than  the 
magnitudes  in  the  18  to  25  gm  F^O  rotational  band. 
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D.  RADIATION  FROM  NOCTILUCENT  CLOUDS 


General  Characteristics  of  Noctilucent  Clouds 

Notcilucen  clouds  have  been  observed  scientifically  since  the  late  19th 
century,  and  their  characteristics  have  slowly  been  elucidated  as  more 
advanced  probing  methods  and  a  larger  number  of  observing  sites  have  been 
utilized  during  the  past  decade.  However,  many  unanswered  questions  re¬ 
main  such  as  the  exact  composition  and  origin  of  the  particles  in  these  clouds, 
the  reasons  for  their  relatively  large  number  density,  and  their  actual  fre¬ 
quency,  duration  and  geographical  extent.  J  he  prevalent  technique  of  visual 
sighting  only  provides  a  partial  source  of  information  about  the  actual  per¬ 
manence  of  noctilucent  clouds  since  the  conditions  for  their  visual  observation 
are  restricted  to  a  relatively  narrow  period  of  time  and  geographical  location 
(Fogle  and  Haurwitz,  1966). 

Ihe  present  state  of  knowledge  of  noctilucent  clouds,  as  related  to  the  problem 
of  infrared  brightness  with  which  we  are  concerned,  can  be  described  as 
follows  (Fogle  and  Haurwitz,  1966;  Fioeco  and  draws,  1969): 

(a)  Average  height:  about  82  km  (observed  range  72  to  92  km) 

(b)  Average  thickness:  about  2  km  (probable  range  0.  5  to  5  km) 

(c)  Latitude  of  observations:  45°  to  80  N  and  S,  most  frequent 
at  about  60°) 

(d)  Seasonal  characteristic:  sharp  summer  maximum 

(e)  Spatial  extent:  10 4  to  more  than  4  x  106  cm2. 

(f)  Duration  of  visual  observation:  several  minutes  to  several  hours. 


12  5 


(g)  Particle  size  distribution  equation:  N  =  Kd  ^ 

(Fogle  and  Ilaurwitz,  1966;  Ilemenway  et  al.  ) 
where:  N  is  the  number  density  of  particles  larger 

_3 

than  d  (typically  N  -  200  cm  for  d  >  0.05  gm) 
d  is  particle  diameter, 
p  varies  between  3  and  4. 

(h)  Particle  size  range:  0.05  gm  to  about  0.5  gm  diameter. 

(i)  Estimated  particle  composition:  probably  meteoric  for 
particles  smaller  than  about  0.  2  or  0.3  gm,  and  mostly 
ice  for  larger  particles  (by  nucleation). 

(j)  Temperature  range  of  mesosphere  between  72  and  92  km: 

135°K  to  more  than  200CK.  Probably  temperature  in  the 
presence  of  NLC:  135  to  160CK. 

(k)  Longitudinal  variations.  Not  clearly  established,  although 
there  is  some  evidence  of  regions  of  preferential  sightings 
such  as  Scandinavia  and  north-western  U.  S.S.R. 

The  relatively  high  particle  number  concentration  (see  above  and  Figure  52) 
and  the  possibility  of  sighting  with  long  slant  paths  (limb  view)  from  outside 
the  atmosphere  indicates  a  possible  source  of  significant  infrared  radiation 
from  these  noctilucent  clouds. 
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2.  Rayleigh  Scattering 

For  the  range  of  particle  sizes  cited  above,  only  the  Rayleigh  type  scattering 
has  to  be  considered.  The  Rayleigh  scattering  equation  as  related  to  the 
geometry  of  this  problem  is: 
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2  2 
9-  Vp 


n.  =  - 4 

X  Xo  2\ 


nL 


m2  -  1 


m“  +  1 


\  2 


(1  +  cos2  0) 


where: 

N  =  scattered  intensity 

X 

N  =  incident  radiation  intensity 

Xo 

n  =  particle  number  concentration 

Vp  =  particle  volume 

\  =  wavelength 

m  =  index  of  refraction  of  particle 

0  =  scattering  angle 

L  =  path  length 


(56) 


The  integral  size  distribution  of  Figure  52  shows  the  number  of  particles 
larger  than  a  given  diameter  d  (Blandari  et  al. ,  1968;  Farlow,  1968). 


If  a  vertical  path  length  of  1  km  is  assumed  and  a  lower  cutoff  diameter  of 
0.  05  /^m  is  used,  N  =*  200  for  d  >  0.05  M m  (d  =^0.5  nm). 

ID  uA 

Starting  from  equation  (55)  and  assuming  the  worst  case  condition  of  p  =  3, 
obtains : 

3 

K  =  Ndp  =  200  (5  x  10"6)  =  2.5  xlO"14  (57) 

,3 

2  11  d 
Since  i  a  V  n  and  V  = 

P  P  6 


one 
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(-pKd'p_1)  =  -2. 06  x  10“14  d2 


5  x  10 
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,,  2  ~  ^  _  nc  in-14  ,2  ,  ,  Q  .  in-28  3 

V  n  =  -  r  2,06  x  10  d  dd  =  8. 5  x  10  cm 

p  J , 


5  x  10 


-5 


2 

Replacing  this  value  of  V  n  in  equation  (56)  one  obtains: 

P 


L  /  m  -  1 

N  =  3. 78  x  10  N  . 

X  xo  x4 


m2  +  1 


(i  +  cos  e) 


(58) 


The  maximum  Rayleigh  scattering  contribution  would  be  obtained  at  \  =  5  m m , 
within  the  spectrum  of  interest.  A  worst  case  calculation  is  based  on  the 
following  values: 

-4 

X  =  5  /nm  =  5  x  10  cm 

m  =  1.  247  for  ice  (Irvine  and  Pollock,  1968) 

9  =  0°  or  180°  (maxima) 

7 

L  =  3  20  km  =  3.  2  x  10  cm  (worst  case  assuming  a  2  km  thickness 

viewed  tangentially  from  outside  the  earth's  atmosphere). 


N  is  determined  by  the  additive  contributions  from  the  solar  and  terrestrial 
Xo 

irradiances  at  5  nm,  estimated  from  theoretical  computations  in  Section  III. 

N  =  7  x  10~4  +  4  x  10“4  watts/cm2  -  Mm1  =  1.  1  x  10  3  watts/cm2-Mm 
\o 

Replacing  these  values  in  equation  (58)  one  obtains: 


,-26 


-3 


3.2  x  10 


3. 78  x  10  xl.lxlO  x 


5  x  10 


-16 


-9  ,2 

1.98  x  10  1  watts /cm  -jum 
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For  one  steradian: 


N 


X 


- 10  2 

=  1.57  x  10  watts/cm  -ster-pm 


The  value  of  for  the  parameters  chosen  above  would  vary  as  shown  in 
Table  3  as  a  function  of  wavelength,  on  the  basis  of  a  composite  Sun-earth 
irradiance  at  about  140  km  above  the  earth's  surface,  assuming  extreme 
magnitudes  of  irradiance  from  the  earth-atmosphere  system. 


Table  3.  Infrared  Scattering  from  INoctilucent  Clouds 
Path  Length  3  20  km  -  Limb  View 


P m 

9 

N  ,  watts/cm“-ster-pm 

A. 

3 

7. 8  x  10“9 

4 

8. 75  x  10"10 

5 

1. 57  x  10“10 

8 

1. 52  x  10-11 

10 

6. 8  x  10"11 

15 

2.  5  x  10  1 

It  thus  appears  that  for  wavelengths  shorter  than  5  jum  the  scattered  intensity 
grows  very  rapidly,  almost  one  order  of  magnitude  for  every  micrometer 
of  decreasing  wavelength. 


The  magnitudes  of  radiance  at  wavelengths  greater  than  5  pm  are  on  the 
-10  -12  2 

order  of  10  to  10  watts /cm  -ster-jum,  respectively.  These  magnitudes 
are  considerably  lower  than  those  expected  from  gaseous  emission  processes 
from  about  5  pm  to  10  pm,  and  are  lower  in  magnitude  than  thermal  emission 
from  interplanetary  dust.  Therefore,  infrared  scattering  from  noctilucent 
clouds  can  be  neglected. 
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3.  Thermal  Emission  from  Noctilucent  Clouds 


The  problem  of  determining  the  thermal  emission  from  noctilucent  clouds  is 
uncertain  since  the  actual  particle  temperature  and  composition  are  unknown. 
The  calculation  of  this  temperature  is  complicated  by  the  fact  that  Rayleigh 
particles  are  cooler  than  the  black  body  temperature  since  their  absorption 
cross  section  is  less  than  their  emission  cross  section,  the  actual  mean 
ratio  of  these  cross  sections  being  dependent  on  the  composition  of  the 
particles.  Even  if  it  is  assumed  that  these  are  ice  particles,  the  problem 
is  not  entirely  resolved  since  at  an  altitude  of  about  80  km  thermal  conduction 
by  air  molecules  still  plays  a  role  (probably  in  the  direction  of  cooling  the 
particles).  In  any  case,  if  the  particles  are  indeed  ice  spheres,  their 
temperature  will  most  probably  be  lower  than  200 °K.  Most  authors  expect  their 
temperature  to  be  as  low  as  130  to  150°I<  (Eogle  and  Haurwitz,  1066). 


The  thermal  emission  from  noctilucent  clouds  can  be  calculated  from  the 
following  integral  (Kaiser's  (1968)  nomenclature  is  being  utilized  in  this 
equation): 


where : 


a^  and  a^  are  the  minimum  and  maximum  particle  radii  of  the 
integral  particle  size  distribution  (see  equation  (55)). 

T3  1 

n(a)  da  =  -K  *  as  defined  before  (scattering  analysis) 


C  (a) 
A. 

where 


=  4/3  rra  6  =  absorption  or  emission  cross  section, 

X 

§  is  the  Rayleigh  coefficient  defined  as  follows: 

X 
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u  v  2  2 

X  (e  +2)  +  fj 

_ 2  —  _ 

where  e  =  n  -  n'  and  n  =  2  n  n',  where  s  is  the  dielectric  constant,  ri  is  the 
resistivity,  n  and  n'  are  the  real  and  the  imaginary  parts  of  the  particle  index 
of  refraction  at  wavelength  B^(T)  is  the  Planck  function  (at  wavelength  x 
and  temperature  T).  L  is  the  path  length  within  the  emitting  cloud. 

Performing  the  above  replacements  we  obtain: 

a 

N  =4/3 rr  §  B  (T)  LK  f  a  "  da  (60) 

K  A  X  P  Ja 

1 

Equation  (60)  was  incorporated  in  a  computer  program  to  calculate  N  as  a 
function  of  \  for  various  dust  size  distribution  and  composition  models,  as 
well  as  a  function  of  the  particle  temperature  T.  Calculations  were  per¬ 
formed  for  two  sets  of  cloud  characteristics.  A  "worst  case"  path  length  of 
320  km  was  assumed  which  represents  a  2  km  thick  cloud  viewed  tangentially 
(limb  view)  from  outside  the  earth's  atmosphere.  Particle  concentrations  and 

size  distribution  were  assumed  to  be  the  same,  200  part.  /cm°  and  a  = 

o 

0.  025/i  to  Uj  =  0.250  Mm,  respectively.  Temperatures  of  135°K  and  160°K 
ind  size  distribution  powers  of  3.  5  and  3.0  (Fogle  and  Haurwitz,  1966)  were 
assumed  as  being  representative  and  upper  limit  cases,  respectively.  Three 
models  of  particle  composition  were  used:  (1)  an  all  dust  cloud  composed  of 
dirty  quartz,  (2)  an  all  ice  sphere  cloud,  and  (3)  a  composite  ice -dust  model 
where  the  small  particles  below  0.  05  gm  were  dust  and  the  iarger  particles 
s  0.05  gm  were  ice  spheres  on  the  basis  of  Ilemenway's  determinations. 

Substitution  of  the  above  physical  parameters  into  Equation  (60)  lead  to  the 
spectra  presented  in  figures  53,  54  and  55.  Figure  53  depicts  the  spectra 
for  the  upper  limit  (1)  and  typical  (2)  cases  of  a  noctilucent  cloud  composed 
exclusively  of  dust,  [he  results  in  figure  54  are  based  on  a  cloud  composed 
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SPECTRAL  RADIANCE  N.  (WATTS/CM^  STER  uM) 


Figure  53.  Variation  in  Spectral  Radiance  for  a  Noctilucent  Cloud 
Composed  Exclusively  of  "Dirty  Quartz".  Curve  1 
Represents  an  Upper  Limit  and  Curve  2  Represents 
the  Typical  Case 
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SPECTRAL  RADIANCE  Nx  (WATTS/CM  STER  -  *iM) 


NOCTILUCENT  CLOUD  EMISSION 
ICE  MODEL  -  LIMB  VIEW 
PATH  LENGTH  320  KM 


WAVELENGTH  X  (yM) 


Figure  54.  Variation  in  Spectral  Radiance  for  a  Noctilucenf  Cloud 

Composed  Exclusively  of  Ice  Particles.  Curve  j  Repre- 
sents  an  Upper  limit  and  Curve  2  Represents  the  Typical 
C  a 


SPECTRAL  RADIANCE  N,  (WATTS/CM 


exclusively  of  ice  particles,  and  Figure  55  shows  the  spectra  for  the  ice-dust 
model  cloud,  a  best  estimate.  The  variations  in  the  magnitudes  of  the  spectra 
with  wavelength  are  due  to  strong  reststrahlen  bands  of  the  materials 
assumed. 

The  ice-dust  model  represents  the  most  probable  magnitudes  of  thermal 
emission  to  be  expected  from  noctilucent  clouds.  In  general,  a  magnitude 

“8  2 

of  about  10  watts/ cm "'-ster-^i  would  be  a  reasonable  spectral  radiance. 

It  is  emphasized  again  that  this  represents  the  worst  case  of  viewing  a 
cloud,  edge -on  with  the  entire  fie  Id -of -view  of  a  radiometer  filled  by  the 
cloud.  Had  the  direction  of  view  been  in  the  vertical  sense,  the  path  length 
would  be  about  3  km  instead  of  320  km  for  the  limb  view.  Under  these  con¬ 
ditions,  the  radiance  magnitudes  would  be  reduced  by  two  orders  of  magnitude, 
to  about  10  watts /cm“-ster-/um. 

In  summary,  thermal  emission  from  noctilucent  clouds  is  of  considerable 
importance  when  considering  the  natural  infrared  radiant  background  in  the 
70  to  90  km  altitude  region.  The  estimate  of  10 "8  watts/ cm2 -ster-ptm  is 
definitely  a  competing  magnitude  of  radiance  as  compared  to  the  magnitudes 
expected  from  gaseous  sources  at  these  altitudes.  Limb  radiances  from 
noctilucent  clouds  would  definitely  be  the  dominant  source  at  1 1  jum  and  in  the 
16  to  25  g m  spectral  range.  The  transient  and  unpredictable  nature  of 
noctilucent  clouds  makes  this  source  extremely  important  when  considering 
variations  in  natural  background  sources. 
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appendix  a 

TABULATED  DATA  FOR  ADOPTED 
ATMOSPHERIC  PROFILES 


As  noted  in  Section  II,  this  appendix  contains  tables  of  the  adopted  tempera¬ 
ture,  water  vapor,  ozone,  and  nitrous  oxide  profiles.  Data  are  tabulated 
only  at  key  heights  which  represent  break  points  in  the  particular  profile. 
Ozone  values  above  3  5  km  and  nitrous  oxide  values  above  40  km  may  be 
obtained  from  analytic  formulas  given  in  the  text.  Each  parameter  profile 
has  a  number  code  associated  with  it.  A  summary  of  these  number  codes  is 
contained  in  Appendix  B. 
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Table  A-l.  Data  for  Adopted  Temperature  Profiles 
(Heights  in  km.  Temperatures  in  °K) 


Code:  1 

Standard 

Atmosphere 

Ht.  T 

2 

15° 

Annual 

Ht.  T 

3 

45° 

Jan. 

Ht.  T 

4 

45° 

July 

Ht.  T 

0 

288 

0 

300 

0 

2  72 

0 

294 

11.  0 

217 

2.3 

286 

3.  0 

2  72 

2.0 

285 

20.  1 

217 

2.  5 

287 

10.0 

220 

6.  0 

261 

32.2 

229 

16.  6 

193 

19.  1 

215 

13.0 

216 

47.4 

271 

22.1 

215 

27.  1 

215 

17.  1 

216 

52. 4 

271 

47.  5 

270 

32.2 

219 

27.  1 

228 

61. 6 

253 

51.  5 

270 

47.4 

266 

32.  2 

238 

80.  0 

181 

59.  7 

254 

52.4 

266 

47.4 

2  76 

90.  0 

181 

82.2 

177 

64.  7 

242 

52.  4 

2  76 

90.  0 

177 

85.  1 

200 

62.6 

165 

90.  0 

200 

82.  0 

165 

90.  0 

165 
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Table  A-l.  Data  for  Adopted  Temperature  Profiles  (Continued) 
(Heights  in  km,  Temperatures  in  °K) 


Code:  5 

6 

7 

8 

75° 

75° 

75° 

75° 

Jan.  (warm) 

Jan.  (cold) 

Jan.  (mean) 

July 

Ht.  T 

lit.  T 

Ht.  T 

Ht.  T 

0 

249 

0 

249 

-  0 

249 

0 

2  78 

1.  5 

2  54 

1.  5 

254 

1.  5 

2  54 

2.  5 

2  72 

8.  9 

215 

8.  5 

215 

8.  5 

215 

9.  5 

226 

11.0 

222 

17.  6 

202 

11.5 

214 

13.5 

230 

17.  1 

222 

2  5.  1 

197 

19.  0 

208 

23.6 

230 

2  5.  1 

226 

30.  2 

197 

30.  1 

208 

30.2 

23  8 

30.  2 

234 

52.  5 

243 

3  7.  5 

212 

40.  0 

260 

52.  5 

263 

90.  0 

210 

3  7.  5 

228 

50.  0 

270 

90.  0 

210 

42.  5 

242 

60.0 

26  5 

47.  5 

2  50 

70.  0 

220 

52.  5 

253 

80.  0 

1  70 

90.  0 

210 

85.  0 

150 

90.  0 

145 
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Table  A-l.  Data  for  Adopted  Temperature  Profiles  (Continued) 
(Heights  in  km,  Temperatures  in  °K) 


Code: 
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Table  A-2.  Data  for  Adopted  Water  Vapor  Mixing  Ratio  Profiles  (Continued) 
(Heights  in  km;  Mixing  Ratios  in  ppm) 
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Table  A-3.  Data  for  Adopted  Ozone  Profiles  (Continued) 

(Heights  in  km;  Ozone  Values  in  ppm  by  Weight) 
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Table  A3.  Data  for  Adopted  Ozone  Profiles  (Continued) 

(Heights  in  km;  Ozone  Values  in  ppm  by  Weight) 
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Table  A-3.  Data  for  Adopted  Ozone  Profiles  (Continued) 

(Heights  in  km;  Ozone  Values  in  ppm  by  Weight) 
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Table  A-4.  Data  for  Adopted  Nitrous  Oxide  Profiles 

(Heights  in  km;  N9O  Values  are  Fractional 
Volume  Abundances) 


Code: 

1 

9 

La 

3 

4 

Strat. 

z2  - 1 

(cm“sec  ) 

104 

102 

103 

105 

Ht. 

N2° 

n2o 

N2° 

N2o 

0 

2.  5  x  10'7 

2. 5  x  10'7 

2. 5  x  10" 7 

2.  5  x  10"7 

10 

2.  5  x  10"7 

2.5x10  7 

2.5x10  7 

2.5x10  7 

15 

2.4  x  10'7 

7.5x10  8 

2. 0  x  10"7 

2.  5  x  10' 7 

20 

2.  1  x  IQ'7 

2. 3  x  10"8 

l. 5  x  10'7 

2.4  x  10' 7 

25 

1. 75  x  10“ 7 

6.  8  x  10" 9 

4.  5  x  10"8 

2.25  x  10'7 

30 

1.25  x  10" 7 

-9 

2,  1-10 

1. 4  x  10"8 

2.  0  x  10'7 

3  5 

3. 8  x  10~8 

6. 1  x  10"10 

4.  1  x  10'9 

1. 6  x  10"7 

40 

1. 1  x  10“8 

1. 9  x  10' 10 

1. 2  x  10"9 

1.15x10  7 
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APPENDIX  B 

CODES  FOR  ATMOSPHERIC  PROFILES 


This  appendix  contains  a  summary  of  the  identification  of  codes  associated 
with  the  tabulated  data  for  adopted  atmospheric  profiles  in  Appendix  A. 
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Table  B-l.  Codes 


T 

II^O  (Lower  Atm) 

O3  CC>2  UgO  (Upper  Atm) 

Standard 

Atmosphere 

1 

1 

111  1 

Annual  Average 
15° 

2 

2 

2 

January  45° 

3 

3 

3 

July  45° 

4 

4 

4 

Winter  (mean) 
75° 

5 

5 

5 

July  75° 

6 

6 

6 
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Table  B-l.  (Coni.) 


H90  (Lower 

Atm.  )  and 

Code 

Annual 

15° 

I  -OW 

1 1 

Annual 

15° 

High 

12 

Summer 

45° 

Low 

13 

Summer 

45 

High 

14 

Winter 

45° 

Low 

15 

Winter 

45° 

High 

16 

Summer 

75° 

-L-pw  i 

17 

Summer 

75° 

High 

\ 

18 

\\  inter 

75° 

i  .ow 

16 

Winter 

75° 

High 

20 

NOTE : 

0)  h9o 

Low  and  High  based  upon  5'r 

and  95"o 

percentile  values. 

Oo  Low  and  1 

f igh  based 

upon 

±2a  values,  or  2.  67.  5T 

percentile. 
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Table  B-l.  (Cont.  ) 


T 

Code 


Winter  (cold  strat.  ) 

75° 

7 

Winter  (warm  strat.  ) 

75° 

O 

o 

Cold  extreme 

9 

Hot  extreme 

10 

H^O  (Upper  Atmos.  ) 


Mean 

Dry 

Wet 


N2° 

Code 
1 
2 

3 

4 

NOTE:  K  refers  to  vertical  exchange  coefficient 
in  stratosphere. 


K  (cm2  sec  *) 
z 

4 

10 


10' 


iol 


10* 


Code 

J. 

2 

3 


162 


APPENDIX  C 

SELECTION  OF  MODEL  ATMOSPHERES  FOR 
R A DIA  NC’  E  C A  LC  U LA T IO N S 

This  appendix  contains  a  selection  of  47  models  for  radiance  calculations  in 
the  lower  atmosphere.  These  model  atmospheres  are  representative  of  the 
spatial  and  temporal  variations,  means  and  extremes  of  atmospheric  condi¬ 
tions  that  are  likely  to  be  encountered  about  the  earth. 

The  identification  of  these  codes  is  presented  in  Appendix  B,  and  these  codes 
are  used  as  input  values  to  the  computer  program  which  compiles  the  model 
atmospheres. 

It  is  noted  that  any  combination  of  codes  presented  in  Appendix  B  can  be  used 
to  develop  a  wide  variety  of  model  atmospheres.  The  47  models  selected  here 
are  only  a  small  but  representative  fraction  of  the  total  number  of  possibilities 
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Table  C  -  1 .  Selection  of  Models  for  Radiance  Calculations 
in  the  Lower  Atmosphere 


Number  of 
Profile 


h2o 

(trop) 


h2o 

(strat) 


EXTREME  TEMPERATURE  PROFILES 
9  9  1  1 

10  10  1  1 


H20  (trop.  )  AND  Og  VARIABILITY 


11 

2 

11 

1 

12 

2 

12 

1 

13 

4 

13 

1 

14 

4 

14 

1 

15 

3 

15 

1 

16 

3 

16 

1 

17 

6 

17 

1 

1 

1 

3 

4 

5 

6 
5 
5 

1 

1 

11 

12 

13 

14 

15 

16 
17 
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Table  C-l.  (Cont.  ) 


Number  of 
Profile 

T 

h2o 

(trop.  ) 

h2o 

(strat) 

°3 

II20  (trop.  )  AND  Og  VARIABILITY  (Cont.  ) 

18 

6 

18 

1 

18 

19 

5 

19 

1 

19 

20 

5 

20 

1 

20 

EXTREME  CONDITIONS,  ALL  PARAMETERS 

21 

10 

12 

3 

20 

22 

10 

12 

3 

12 

23 

9 

19 

2 

17 

CLIMATOLOGICAL  MEANS  WITH  MOIST  STRATOSPHERE 

24 

1 

1 

3 

1 

25 

2 

2 

3 

2 

26 

3 

3 

3 

3 

27 

4 

4 

3 

4 

28 

5 

5 

3 

5 

29 

6 

6 

3 

6 

30 

7 

5 

3 

5 

31 

8 

5 

3 

5 
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Table  C-l.  (Cont.  ) 


Number  of 
Profile 

T 

h2o 

(trop.  ) 

h2o 

(strat) 

°3 

CLIMATOLOGICAL  MEAN  WITH  DRY  STRATOSPHERE 

32 

1 

1 

2 

1 

33 

2 

2 

2 

2 

34 

3 

3 

2 

3 

35 

4 

4 

2 

4 

36 

5 

5 

2 

5 

37 

6 

6 

2 

6 

38 

7 

5 

2 

5 

39 

8 

5 

2 

5 

STANDARD  ATMOSPHERE  TEMPERATURE  PROFILE 
EXTREMES  OF  H20  AND  Og 


40 

1 

12 

3 

12 

41 

1 

12 

3 

20 

42 

1 

19 

2 

17 

43 

1 

1 

1 

12 

44 

1 

1 

1 

20 

45 

1 

1 

1 

17 

46 

1 

19 

2 

1 

47 

1 

12 

3 

1 
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APPENDIX  D 

ATMOSPHERIC  MODEL  FOR  NITRIC  ACID 


/  JUSTIFICATION  FOR  THE  MODEL 


Recent  measurements  made  during  balloon  flights  flown  by  the  University  of 
Denver  have  shown  the  existence  of  nitric  acid  in  the  earth's  atmosphere  (Mur- 
cray,  et  al,  1968).  Identification  of  nitric  acid  was  made  using  a  balloon-borne 
spectrograph  at  an  altitude  of  30  km  to  observe  the  solar  spectrum  at  zenith 
angles  near  90°.  Identification  was  made  by  the  absorption  anomaly  at  7.  5  microns. 

Nitric  acid  also  has  a  strong  absorption  band  near  11.2  microns,  which  is 
important  to  the  present  atmospheric  studies.  The  present  investigation  has 
been  to  assess  the  radiation  from  the  limb  from  5  to  25  microns  with  emphasis 
on  the  window  regions  where  the  horizon  background  is  low  for  enhancement 
of  detection  of  cold  objects  in  space.  The  previously  used  atmospheric  models 
have  consistently  shown  a  transparent  atmospheric  window  in  the  spectral 
region  from  10  to  12  microns.  This  discrepancy  between  the  previously  exist¬ 
ing  models  and  recent  measurements  has  made  it  necessary  for  the  inclusion 
of  a  preliminary  model  for  nitric  acid  absorption;  even  a  crude  model  to  give 
first-order  effects  is  sufficient  for  the  design  of  a  space  experiment  and 
should  be  more  accurate  than  its  omission  from  the  limb  calculations. 

A  subsequent  balloon  flight  was  made  by  the  University  of  Denver  with  a 
spectral-radiometer  system  capable  of  measuring  the  absolute  radiance  of 
the  earth  and  its  atmosphere  in  the  wavelength  region  from  5  to  20  microns 
at  selected  elevation  angles.  Nine  spectra  have  been  reported  for  the  spectral 
interval  from  10  to  13  microns  for  three  different  altitudes  and  for  three  dif¬ 
ferent  elevation  angles  near  the  horizon  (Murcray  1969).  An  analysis  of  the 
nine  spectra  was  made  to  give  a  first-order  approximation  for  the  limb 
radiance  calculations. 
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B.  DETERMINATION  OF  THE  MODEL 


Polyatomic  molecules  in  the  gaseous  phase,  such  as  atmospheric  water  vapor, 
ozone  and  methane  exhibit  vibration -rotation  spectra  at  infrared  frequencies. 
It  has  been  shown  that  for  pressures  and  temperatures  encountered  in  the 
lower  atmosphere,  a  good  approximation  to  the  absorption  can  be  obtained 
from  the  statistical  absorption  model  (Goody,  1952),  given  by: 


t  =  exp 


'  -  3x 

./T+aE 


(Dl) 


where 


t  =  transmission 

o  =  2tt  AyL 

p  v 

S  w 

X  2rr  Av  L 

S  =  average  strength  of  the  absorption  lines 

v  =  average  frequency  spacing  between  lines 

w  =  optical  path  of  absorber  reduced  to  standard  conditions 


(D2) 

(D3) 


average  width  of  absorption  lines,  measured  as  half -width 
at  half  peak  intensity.  The  half-width  has  been  shown  to  be 
directly  proportional  to  pressure  and  inversely  proportional 
to  the  square  root  of  temperature.  For  the  atmospheric 
paths  involved  in  the  Denver  measurements,  the  temperature 
variation  was  small,  therefore  the  temperature  contribution 
to  half -width  can  be  neglected. 


For  long  slant  paths  through  the  atmosphere,  where  the  line  width  is  very 
small  compared  to  the  line  strength,  then  absorption  by  an  increase  in 
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optical  patli  will  predominently  result  from  the  wings  of  the  rotational  lines 
and  the  statistical  model  can  be  approximated  by  the  following  equation: 


t  =  exp 


(1/2  82  x) 


(1)4) 


The  optical  slant  path  through  the  atmosphere  can  be  approximated  by  a 
C  urtis -Godson  slant  path  integration  (Curtis,  1952).  For  a  uniformly 
distributed  molecule,  the  optical  path  should  be  proportional  to  the  density 
times  pressure,  integrated  over  the  path,  and  transmission  should  be 
approximately  given  by: 


t  =  exp 


K  /J  o  p  d  L 


where  p  is  the  atmospheric  density 


(1)5  i 


1'  or  the  altitudes  represented  by  the  Denver  measurements,  the  atmosphere  li¬ 
nearly  in  local  thermodynamic  equilibrium  and  the  net  radiation  received  at 
the  detector  is  given  by: 

N.  =  J  B.  d  t  (1)6) 


where 

N  is  the  net  atmospheric  radiance 

A. 

B  is  the  Planck  radiation  function  at  the  temper atui  e  of  the 
atmosphere  and  wavelength  measured 

Atmospheric  transmission  was  estimated  at  peak  absorption  wavelength  of 
11.2  microns  for  each  of  the  nine  University  of  Denver  measurements  by 
dividing  measured  radiance  by  the  Planck  function  calculated  for  the  average 
temperature  in  the  line  of  sight. 
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figure  Dl  shows  a  plot  of  the  quantity 

~  loge  * 

V  '  w  p  d  L 

vs  the  effective  altitude  of  the  observation. 

The  effective  altitude  was  taken  as  2.  8  km  above  the  minimum  altitude  of  the 
ray  path.  If  nitric  acid  were  uniformly  distributed  over  the  altitudes  for  the 
measured  data  then  the  above  parameter  would  be  a  constant,  independent  of 
altitude.  However,  Figure  Dl  indicates  that  the  mixing  ratio  is  increasing 
with  altitude.  The  figure  shows  a  nearly  linear  increase  in  the  parameter 
which  requires  that  mixing  ratio  be  proportional  to  the  square  of  the  altitude. 
An  iterative  inversion  technique  was  attempted  to  fit  the  following  expression 

n  T  3  / 2  o 

KjVpE)  (^2)  //  d  L 

o 

where  Pq  and  are  taken  as  normal  temperature  and  pressure  and  Z,  is 
the  altitude  above  the  earth. 

The  result  of  this  inversion  gave  a  value  of  K  =  0.  84  x  10  ~12,  which  is  inter¬ 
preted  as  a  product  of  the  nitric  acid  absorption  parameter  times  a  mixing 
ratio  parameter. 

C.  CONCLUSIONS  AND  LIMITATIONS  OF  THE  MODEL 

The  radiance  values  calculated  from  the  model  were  compared  with  the 
measured  values  at  11.2  microns.  Figure  D2  shows  the  comparison  for 
each  of  the  nine  measurements,  and  shows  fair  agreement  except  for  one 
value.  This  one  extreme  value  was  obtained  for  the  lowest  measurement 
with  the  instrument  looking  down  into  the  'dmosphere  at  the  zenith  angle 
of  93.5  .  ri  his  measurement  represents  a  minimum  ray  path  altitude  of 


1/2 


(D7) 


t  =  exp 


{ 
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approximately  6.4  km.  This  extreme  value  leads  to  the  conclusion  that  the 
nitric  acid  mixing  ratio  falls  of  much  more  rapidly  than  that  given  by  the 
model  for  altitudes  below  10  km. 

If  the  atmospheric  absorption  and  emission  does  sensibly  behave  like  the 
model  assumed  above  and  if  the  line  strength  is  reasonbly  independent  of 
pressure,  then  the  relative  density  of  nitric  acid  in  the  atmosphere  will 
be  given  as: 


°IIN()3  K'  Z  pAir  (,) 

where 

"  H.\On  =  relative  distribution  of  nitric  acid 
0 

,  .  =  density  of  the  air 

Air  J 

V.  -  altitude 

K'  =  normalization  constant  chosen  to  give  a  relative  distribution 
of  one  at  the  altitude  of  peak  nitric  acid  density. 

Figure  L)3  shows  the  relative  distribution  for  the  nitric  acid  model.  The  pea 
density  occurs  at  an  altitude  of  approximately  12  km. 

Dr.  Williams  of  the  University  of  Kansas  has  considered  the  results  of  this 
inversion  ana  has  suggested  the  following  explanation  (Williams,  1 D 6 D ) . 

Nitric  acid  may  be  formed  from  ozone  in  the  atmosphere.  The  peak  of  the 
ozone  distribution  has  been  observed  at  altitudes  of  approximately  20  to  30 
kilometers.  Doing  a  heavy  molecule,  nitric  acid  could  diffuse  to  lower  alti¬ 
tudes  until  it  came  into  contact  with  liquid  water,  whose  propensity  for 
adsorption  of  nitric  acid  is  high.  Water  droplets  are  not  likely  to  be  en¬ 
countered  above  the  tropopause.  The  atmospheric  tropopause  at  mid¬ 
latitudes  is  at  about  10  to  12  kilometers. 
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The  above  explanation  indicates  that  the  model  determined  for  nitric  acid 
may  be  realistic  and  appears  to  be  adequate  for  a  first-order  approximation 
f°r  the  design  of  an  experiment.  However,  it  must  be  realized  that  the 
model  was  synthesized  from  a  minimum  of  experimental  evidence  and  that 
the  assumption  of  a  different  absorption  model  may  give  an  entirely  different 
distribution;  also,  the  model  has  been  extrapolated  to  high  altitudes  and 
rapidly  loses  validity  above  30  to  3  5  km. 
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Figure  D1 .  Test  of  Hypothesis  of  I'niform 
Distribution  of  UNO,, 
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Figure  D2.  Comparison  between  Measurements 
and  Model  At  11.2  Microns 
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* 


Figure  D3.  Atmospheric  Distribution  of  Nitric  Acid 


t 
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APPENDIX  L 


SPECTRAL  RADIANCE  FOR  47  SELECTED 
LOWER  ATMOSPHERIC  MODELS 


This  appendix  contains  spectra  calculated  for  the  47  atmospheric  profiles 
tabulated  in  Appendix  C,  Table  C-l.  For  each  profile,  the  profile  number 
is  given  together  with  the  profile  code  explained  in  Appendix  C.  For  each 
profile,  spectra  are  given  for  scale  heights  from  10  to  80  kilometers,  in 
increments  of  10  kilometers.  The  scale  height  for  each  spectrum  is 
labeled.  The  abscissa  is  given  as  spectral  wavelength  in  microns  and  the 

_  9  _  i  _  ] 

ordinate  is  given  as  spectral  radiance  in  watts  cm  “  steradian  micron 
and  is  plotted  on  seven  cycle  semilogarithmic  paper. 

Figures  Ell  to  E39  do  not  give  radiance  beyond  10.  5  microns  because  the 
same  information  is  contained  in  other  spectra.  For  each  model,  reference 
is  made  to  the  appropriate  model  for  longer  wavelengths. 


176 


Profile  No.  4  (4-4-1-4) 
Horizon  Spectral  Radiance 
Mid-latitude  Summer  Atmosphere 


Figure  E5 

Profile  No.  5  (5-5-1-5) 
[orizon  Spectral  Radiance 
Artie  Winter  Atmosphere 


10  Km 


Figure  E6 

Profile  No.  6  (6-6-1-6) 
Horizon  Spectral  Radiance 
Arctic  Sumner  Atmosphere 
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Profile  No.  8  (8-5-1-5) 

Horizon  Spectral  Pvadiance 
Arctic  Atmosphere  -  Warm  Stratosphere 


Spectral  Radiance 
Cold  Atmosphere 


Figure  Ell 

Atmospheric  Model  No.  11  (2-11-2-11) 
Horizon  Spectral  Radiance 
Tropical-Dry  Troposphere 

(From  10.5  to  25  microns,  see  Model  No.  2) 
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6  8  10  12 
WAVELENGTH  (MICRONS) 


Figure  E12 

Atmospheric  Model  No.  12  (2-12-1-12) 
Horizon  Spectral  Radiance 
Tropical-Moist  Troposphere 

(From  10.5  to  25  microns,  see  Model  No.  2) 
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RADIANCE  (W/CM  -SR-u> 


Figure  E13 


Atmospheric  Model  No.  13  (4-13-1-13) 
Horizon  Spectral  Radiance 
Mid-latitude,  Summer,  Dry  Troposphere 


(From  10.5  to  25  microns,  see  Model  No.  4) 
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6  8  10  1 
WAVELENGTH  (MICRONS) 


Figure  E14 

Atmospheric  Model  No.  14  (4-14-1-14) 

Horizon  Spectral  Radiance 
Mid-latitude,  Summer,  Moist  Troposphere 


(From  10.5  to  25  microns,  see  Model  No.  4) 


6  8  10  12 
WAVELENGTH  (MICRONS) 


Figure  E15 

Atmospheric  Model  No.  15  (3-15-1-15) 
Horizon  Spectral  Radiance 
Mid-latitude,  Winter,  Dry  Troposphere 

(From  10.5  to  25  microns,  see  Model  No.  3) 
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6  8  10  12 
WAVELENGTH  (MICRONS) 


Figure  E16 

Atmospheric  Model  No.  16  (3-16-1-16) 
Horizon  Spectral  Radiance 
Mid-latitude ,  Winter,  Moist  Troposphere 

(From  10.5  to  25  microns,  see  Model  No.  3) 
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Figure  E17 

Atmospheric  Model  No.  17  (6-17-1-17) 
Horizon  Spectral  Radiance 
Arctic  Summer,  Dry  Troposphere 

(From  10.5  to  25  microns,  see  Model  No.  6) 
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Figure  E18 

Atmospheric  Model  No.  18  (6-18-1-18) 
Horizon  Spectral  Radiance 
Arctic  Summer  -  Moist  Troposphere 

(From  10.5  to  25  microns,  see  Model  No.  6) 
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6  8  10  12 
WAVELENGTH  (MICRONS) 


Figure  E19 

Atmospheric  Model  No.  19  (5-19-1-19) 
Horizon  Spectral  Radiance 
Arctic  Winter  -  Dry  Troposphere 

(From  10.5  to  25  microns,  see  Model  No,  5) 
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RADIANCE  (W/CWT-SR-n) 


Figure  E20 

Atmospheric  Model  No.  20  (5-20-1-20) 
Horizon  Spectral  Radiance 
Arctic  Winter  -  Moist  Troposphere 

(From  10.5  to  25  microns,  see  Model  No.  5) 
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Figure  E21 

Atmospheric  Model  No,  21  (10-12-3-20) 
Horizon  Spectral  Radiance 
Extreme  Conditions  -  Maximum  Radiance 

(From  10.5  to  25  microns,  see  Model  No.  10) 
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RADIANCE  (W/CM  -SR-4 ) 


6  8  10  12 


WAVELENGTH  (MICRONS) 
Figure  E22 

Atmospheric  Model  No.  22  (10-12-3-12) 
Horizon  Spectral  Radiance 
Extreme  Conditions  -  Maximum  Radiance 

(From  10.5  to  25  microns,  see  Model  No.  10) 
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6  8  10  12 


WAVELENGTH  (MICRONS) 

Figure  E23 

Atmospheric  Model  No.  23  (9-19-2-17) 
Horizon  Spectral  Radiance 
Extreme  Conditions  -  Minimum  Radiance 

(From  10.5  to  25  microns,  see  Model  No.  9) 
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RADIANCE 


SE 


Figure  E25 

Atmospheric  Model  No.  25  (2-2-3-2) 
Horizon  Spectral  Radiance 
Tropical,  Moist  Stratosphere 

(From  9  to  25  microns,  see  Model  No.  2) 
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Figure  E26 

Atmospheric  Model  No.  26  (3-3-3-3) 
Horizon  Spectral  Radiance 
Mid-latitude  Winter,  Moist  Stratosphere 

(From  9  to  25  microns,  see  Model  No.  3) 
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Figure  E28 

Atmospheric  Model  No.  28  (5-5-3-S) 
Horizon  Spectral  Radiance 
Arctic  Winter,  Moist  Stratosphere 

(From  9  to  25  micr  ms ,  see  Model  No.  5) 
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Figure  E29 

Atmospheric  Model  No.  29  (6-6-3-6) 
Horizon  Spectral  Radiance 
Arctic  Summer,  Moist  Stratosphere 


(From  9  to  25  microns,  see  Model  No.  6) 
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5  7  9  11 


WAVELENGTH  (MICRONS) 


Figure  E30 

Atmospheric  Model  No.  30  (7-5-3-5) 
Horizon  Spectral  Radiance 
Arctic,  Cold  Moist  Stratosphere 

(From  9  to  25  microns,  see  Model  No.  7) 
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Figure  E31 


Atmospheric  Model  No.  31  (8-5-3 -5) 
Horizon  Spectral  Radiance 
Arctic,  Warm  Moist  Stratosphere 


(From  9  to  25  microns,  see  Model  No.  8) 
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Figure  E32 

Atmospheric  Model  No.  32  (1-1-2-1) 
Horizon  Spectral  Radiance 
U,  S.  Standard,  Dry  Stratosphere 

(From  9  to  25  microns,  see  Model  No.  1) 
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RADIANCE  (W/CM  -SR-w ) 


Figure  E33 

Atmospheric  Model  No.  33  (2-2-2-2) 
Horizon  Spectral  Radiance 
Tropical,  Dry  Stratosphere 

(From  9  to  25  microns,  see  Model  No.  2) 
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Figure  E34 

Atmospheric  Model  No.  34  (3-3-2-3) 
Horizon  Spectral  Radiance 
Mid-latitude  Winter,  Dry  Stratosphere 

(From  9  to  25  microns,  see  Model  No.  3) 
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Figure  E35 

Atmospheric  Model  No.  35  (4-4-2-4) 
Horizon  Spectral  Radiance 
Mid-latitude  Summer,  Dry  Stratosphere 

(From  9  to  25  microns,  see  Model  No.  4) 
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Figure  E36 


Atmospheric  Model  No.  36  (5-5-2-5) 
Horizon  Spectral  Radiance 
Arctic  Winter,  Dry  Stratosphere 


(From  9  to  25  microns,  see  Model  No.  5) 
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Figure  E38 

Atmospheric  Model  No.  38  (7 -5-2 -5 ) 
Horizon  Spectral  Radiance 
Arctic,  Cold  Dry  Stratosphere 

(From  9  to  25  microns,  see  Model  No.  7) 
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Atmospheric  Model  No.  40  (1-12-3-12) 
Horizon  Spectral  Radiance 
U.  S.  Standard  Atmosphere 


Atmospheric  Model  No.  41  (1-12-3-20) 
Horizon  Spectral  Radiance 
U.  S.  Standard  Atmosphere 


Atmospheric  Model  No.  42  (1-19-2-17) 
Horizon  Spectral  Radiance 
U.  S.  Standard  Atmosphere 


on  Spectral  Radianc 
Standard  Atmosphere 


zon  Spectral  Radiance 
Standard  Atmosphere 
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Atmospheric  Model  No.  46  (1-19-2-1) 
Horizon  Spectral  Radiance 
U.  S.  Standard  Atmosphere 


1U  Km 
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Atmospheric  Model  No.  <47  (1-12-3-1) 
Horizon  Spectral  Radiance 
U.  S.  Standard  Atmosphere 


APPENDIX  F 


SPECTRAL  RADIANCES  FOR 
UPPER  ATMOSPHERE  MODELS 


The  final  goal  of  this  study  of  upper  atmosphere  infrared  radiance  is  the 
prediction  of  spectral  radiance  values  at  wavelengths  between  5  and  25  gm 
for  all  species  and  for  both  exoatmospheric  and  endoatmospheric  viewing. 
Figure  FI  shows  spectral  limb  radiances  between  5  and  2  5  g m  for  tangent 
heights  of  60  km,  80  km  and  100  km.  These  radiances  were  computed  by  a 
program  which  is  included  as  a  subroutine  in  the  final  radiance  model  com¬ 
puter  program. 

For  all  contributions  to  the  spectra  from  molecular  vibration  rotation  bands, 
the  total  band  radiances  are  spread  out  with  the  assumption  that  the  rotational 
levels  are  in  thermal  equilibrium  at  a  mean  temperature  computed  for  the 
line  of  sight  determined  by  the  populations  of  the  lower  vibrational  level. 
Using  this  mean  temperature,  rotational  populations  are  computed  and  wave¬ 
lengths  and  relative  intensities  in  emission  of  each  line  are  computed.  The 
effects  of  self-absorption  for  each  line  are  determined  from  the  optical  thick¬ 
ness  through  Eq.  41.  The  sum  of  relative  radiances  corrected  for  optical 
thickness  is  normalized  to  the  total  band  radiance  and  absolute  line  radiances 
computed.  A  further  subroutine  then  computes  the  contribution  of  each  line 
to  the  spectrum. 

In  the  case  of  linear  molecules,  line  positions  are  computed  for  100  lines  in 
each  of  the  P,  Q  and  R  branches,  neglecting  the  second  rotational  constant 
D^.,  and  any  higher  rotational  constants.  The  effects  of  spin  and  l  doubling 
are  neglected.  In  all  cases,  computed  positions  of  lines  are  within  0.  1  gm 
of  accurate  positions.  For  the  water  vapor  6.3  gm  band,  the  positions  and 
line  lengths  of  about  450  of  the  strongest  lines  tabulated  by  Benedict  and 
Calfee  (1967)  are  used.  For  the  ozone  9.  6  gm  band,  mean  line  strengths  at 
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5  cm  *  determined  from  the  tabulation  of  Clough  and  Kneizys  (1965)  are 
used,  rather  than  attempting  to  include  individual  lines.  The  contribution 
to  the  spectra  from  water  vapor  rotational  lines  is  determined  from  the 
dipole  matrix  elements  and  rotational  level  energies  of  about  250  lines 
between  10  and  25  gm.  Data  for  wavelengths  longer  than  12  g m  were 
determined  from  the  Tables  of  Gates,  et  al.  (1964).  Those  between  10 
and  12  gm  were  determined  from  values  of  absorption  coefficients  supplied 
by  Dr.  Ft.  McClatchey  of  AFCFtL. 

The  present  program  computes  spectra  for  a  triangular  slit  function  with  a 
width  at  half  maximum  of  0.25  gm.  The  appropriate  subroutine  may  be  re¬ 
written  for  any  desired  slit  function.  The  approximations  made  in  the  current 
computational  model  restrict  computer  resolution  to  about  0.  1  gm.  At  the 
expense  of  increased  computation  time  and  complexity,  better  resolution  may 
be  obtained. 
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The  infrared  radiance  arising  from  the  earth's  limb  for  altitudes  from  sea  level  to 
500  kilometers  and  within  the  spectral  limits  of  5  to  25  microns  was  determined  from 
theoretical  calculations.  The  atmosphere  was  divided  into  two  altitude  regimes:  the 
lower  atmosphere  (below  70  kilometers)  and  the  upper  atmosphere  (above  70  kilo¬ 
meters).  The  lower  limb  models  were  developed  from  the  University  of  Michigan 
computer  program  with  modifications,  including  the  addition  of  a  model  for  nitric 
acid.  Spectral  radiances  were  determined  as  a  function  of  geometrical,  geographical, 
and  temporal  parameters.  Emphasis  was  placed  on  developing  limb  radiance  models 
for  the  upper  atmosphere,  including  identification  of  active  molecular  species  in 
the  infrared,  the  specification  of  their  abundances, and  their  wavelengths  of  emission. 
In  addition,  particulate  material  sources,  including  infrared  radiation  arising  from 
dust  particles  and  noctilucent  clouds,  and  uncertainties  in  atmospheric  parameters, 
were  considered. 
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